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1.0  INTRODUCTION 


An  explosive  or  propellant,  usually  confined  as  in  a  warhead,  undergoes 
a  chemical  reaction  when  exposed  to  heat.  The  rate  of  this  reaction  is  tem¬ 
perature-dependent  .  The  thermal  response  can  be  expressed  mathematically  by: 


-k 


V2T  + 


c 

p  at 


p  Q  n  e 


-E/RT 


(1-1) 


where  the  first  term  on  the  left  is  the  heat  conduction  term,  the  second  term 
is  the  sensible  energy  accumulated,  and  the  term  on  the  right  is  a  heat  gener¬ 
ation  term.  The  symbols,  k,  p,  C^,  t,  and  T,  carry  their  usual  connotation  of 
thermal  conductivity,  density,  specific  heat,  time,  and  temperature,  respec¬ 
tively.  The  heat  generation  term  also  depends  upon  the  heat  of  decomposition 
Q,  the  collision  number  n,  the  activation  energy  E,  and  the  gas  constant  R. 

As  the  explosive  heats  up,  the  heat  generated  continues  to  increase  until  det¬ 
onation.  Indeed,  for  a  particular  geometry,  if  a  "layer"  of  explosive  reaches 

a  temperature  above  its  critical  temperature  T ,  even  if  all  other  heat 

c 

sources  (such  as  external  fires)  are  suddenly  removed,  the  explosive  will 
still  detonate  as  the  heat  generated  in  the  explosive  cannot  be  dissipated  by 
conduction  through  the  shell.  Instead,  the  heat  feeds  itself,  releasing  more 

heat,  causing  the  reaction  to  release  more  heat  at  a  faster  rate  .  until 

detonation. 

Studies  by  the  Army  have  shown,  for  Army  munitions,  that  propellant 
fires,  though  extremely  hot,  are  not  likely  to  be  a  problem  because  of  their 
short  duration  (less  than  10  seconds).  However,  longer  duration  fires,  such 
as  those  resulting  from  burning  explosive,  sequential  initiation  of  different 
propellant  charges,  or  pool-type  fires  from  fuels  (such  as  JP-4  or  JP-5)  will 
lead  to  cook-off,  as  it  is  commonly  called. 


The  United  States  Navy  first  became  interested  in  the  cook-off  problem 
as  a  result  of  two  disasterous  fires  aboard  the  aircraft  carriers  USS  FOR- 
RESTAL  (July,  1967)  and  USS  ENTERPRISE  (January,  1969).  Both  were  caused  by 
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accidents  on  tlie  flight  deck.  Spilled  aviation  fnel  caused  the  fire  to  spread 
rapidly.  Weapons  lying  on  the  deck  and  hung  under  the  vings  of  nearby  air¬ 
craft  soon  became  enveloped  in  these  fires.  Within  one  to  three  minutes  after 
the  flames  first  became  evident,  the  weapons  began  to  detonate  or  launch  pro¬ 
pul  sively.  The  catastrophic  fire  aboard  the  USS  FORRES TAL  resulted  in  72 
million  dollars  of  damage  to  the  ship,  and  the  deaths  of  133  men.  The  fire 
aboard  the  USS  ENTERPRISE,  while  not  as  catastrophic,  also  resulted  in  loss  of 
life,  injury,  and  ship  damage  [1]. 

Another  carrier  incident  involving  weapons  and  fire  occurred  aboard  the 
USS  ORISKANY.  This  fire,  attributed  to  the  lack  of  environmental  sensing 
fuzes,  was  not  propagated  by  aircraft  fuel  but  did  result  in  cook-off  of 
numerous  flares  causing  extensive  ship  damage  and  loss  of  life.  In  September 
1970,  at  Danang  Air  Base,  an  accidental  fire  erupted  under  an  F-4B  during  air¬ 
craft  refueling  operations.  As  the  fire  was  being  brought  under  control,  a 
Sparrow  warhead  cooked-off,  killing  3  and  injuring  13.  In  1973,  there  were 
two  railroad  boxcar  fires,  one  in  a  railroad  yard  at  Roseville,  California, 
involving  7000  bombs,  and  another  in  a  remote  location  in  Benson,  Arizona.  In 
both  fires,  bombs  detonated  causing  extensive  property  damage,  and  in  the 
Roseville  fire,  injury  also  resulted.  In  November  1975,  a  collision  between 
the  cruiser  USS  BELKNAP  and  the  aircraft  carrier  USS  JOHN  F.  KENNEDY 
occurred.  Ruptured  fuel  lines  on  the  KENNEDY  spewed  fuel  down  the  stacks  and 
across  the  deck  of  the  BELKNAP.  The  resulting  fire  caused  ordnance  cook-off 
[1],  The  most  recent  cook-off  occurred  aboard  the  aircraft  carrier  USS  NIMITZ 
in  the  spring  of  1981  which  resulted  in  the  immediate  deaths  of  14  men  and 
over  60  million  dollars  of  damage. 

In  1969,  the  Weapons  Cook-Off  Improvement  Program  was  established  to 
accomplish  two  specific  objectives:  delay  the  ordnance  reaction  for  at  least 
five  minutes  and  limit  the  extent  of  that  reaction  to  a  deflagration  (case 
rupture  and  burning  of  explosive).  Since  then,  the  five-minute  delay  has  been 
successfully  achieved  primarily  through  the  use  of  ablative  and  intumescent 
coatings  on  the  exterior  skin  of  various  bombs,  rocket  launchers,  mines. 
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missiles,  etc.  Intnmescent  coatings  are  the  material  of  choice  whenever  high 
insolation  efficiency  from  a  limited  film  thickness  is  desired  [2]. 

Many  thermal  materials  have  been  evaluated  and  several  high  thermal  ef¬ 
ficiency  intnmescent  coatings  have  been  found  to  be  viable  and  are  currently 
used  to  protect  rocket  launchers,  missiles,  and  warheads.  While  chemists  have 
become  quite  adept  at  developing  intnmescent  coatings,  this  is  a  result  of 
some  quantitative  physical  and  chemical  understanding,  a  larger  amount  of 
qualitative  physical  and  chemical  understanding,  experience,  trial  and  error, 
and  "black  magic."  What  is  desired  is  that  a  more  quantitative  understanding 
of  the  intumescence  process  be  obtained  by  analyzing  typical  intnmescent  reac¬ 
tion  mechanisms  and  developing  an  appropriate  model  which  predicts  heat  trans¬ 
fer  and  accounts  for  fundamental  chemical,  thermodynamic,  and  physical  pro¬ 
cesses  such  as  kinetics  and  swelling. 
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2.0  THE  PHYSICAL  PROCESSES  AND  CHEMISTRY  OF  INTUMESCENCE 

la  order  for  intumescence  to  occur,  several  distinct  reactions  must 
occur  nearly  simultaneous,  but  in  the  proper  sequence.  Generally,  the  compo¬ 
sition  of  an  intumescent  system  consists  of  the  following  functional  entities: 

(1)  An  inorganic  acid  or  a  material  yielding  acid  at  temperatures 
between  100°  and  250°C, 

(2)  A  polyhydric  material  rich  in  carbon, 

(3)  An  organic  amine,  amide,  or  azo  compound, 

(4)  A  halogenated  material  or  other  suitable  material  which  can  decom¬ 
pose  and  serve  as  a  blowing  agent, 

(5)  A  hydrated  salt  (e.g.,  borax,  sodium  silicate,  aluminum  hydroxide, 
ammonia  phosphate,  etc.),  and 

( 6 )  A  binder. 

Most  formulations  contain  an  example  of  each,  although  in  some  cases,  two  or 
more  functional  groupings  are  included  in  a  single  component. 

The  state-of-the-art  to  date  in  intumescent  systems  includes  several 
basic  categories  of  material  which  can  be  summarized  as  follows: 

(1)  Elastomer-based,  filled  film  which  contains  such  intumescent  solids 
as  hydrates  of  boron  salts,  alkali  silicates,  and  mineral  hydrates 
such  as  perlite  and  vermiculite.  Elastomeric  binders  include  poly¬ 
urethanes,  silicones,  chloroprenes,  polysulfide,  and  mixtures  of 
these  with  rigid  plastics  such  as  phenolics  and  acetals. 
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(2)  Acetal  binders  such  as  melamine,  phenol  formaldehyde,  nrea  formal¬ 
dehyde  filled  with  hydrated  salts,  and  intnmescent  minerals  as 
-above. 

(3)  Anal ine-bi sulfite  addition  products  which  interne see  to  a  char  upon 
heating/burning. 

(4)  Inorganic  coatings  filled  with  intumescent  compounds  including 
hydrated  salts,  minerals,  and  organic  agents  which  oxidize  to  a 
char  such  as  polysachrides  (sugars,  starches)  and  certain  pheno- 
lics. 

(5)  Intumescent  film  forming  materials  such  as  alkali  metal  silicate 
solutions  with  and  without  fillers. 

(6)  Thermoplastic  material  filled  with  gas  producing  compounds  which 
release  vaporous  gases  such  as  freons,  carbon  dioxide,  and  nitrogen 
during  heating  and  produee  an  expanding  mass  of  flame  inhibited 
plastic  until  the  expansion  is  complete  and  the  barrier  gases 
escape.  The  material  may  then  melt,  burn,  or  char,  depending  upon 
composition. 

The  fire  barrier  qualities  of  each  of  the  above  must  be  weighed  individually 
but,  in  general,  the  mechanism  is  quite  similar. 

A  typical  intumescing  coating,  as  it  is  exposed  to  heat,  has  the  fol¬ 
lowing  history  [3].  First,  the  acid  salt*  must  decompose  to  yield  the  dehy¬ 
drating  acid.  This  inorganic  acid  is  capable  of  esterification  of  polyols. 


*  The  intnmescent  system  investigated  in  this  program  is  discussed  in  detail 
in  Section  4.0.  Because  it  contains  a  hydrated  salt,  the  intumescence 
occurs  by  a  different  reaction  mechanism  than  an  acid  salt. 
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Hitt  it,  the  resultant  tcid  removes  the  water  from  tlie  carbonaceous  polyol  at 
temperatures  slightly  above  the  acid  release  temperature.  Amides  (or  amines) 
probably  catalyze  this  reaction  resulting  in  its  completion  at  temperatures 
lower  than  would  be  observed  in  their  absence.  The  acid  acts  as  a  dehydration 
agent  leaving  a  carbonaceous  residue. 

The  "timing"  (that  is,  the  temperature  sequence  of  the  various  compo¬ 
nents)  is  very  important.  The  acid-carbonaceous  polyol  must  melt  prior  to  or 
during  esterification  so  that  when  the  polyol  decomposes  via  hydration  (form¬ 
ing  a  carbon-inorganic  residue),  the  released  blowing  agent  and  the  evolution 
of  other  nonflammable  gases  cause  the  carbonizing  mass  to  puff  and  foam.  One 
of  the  purposes  of  the  organic  amine  or  amide  (such  as  urea,  melamine,  dicyan- 
diamide,  urea  formaldehyde,  etc.)  is  the  release  of  nonflammable  gases  (CO^, 
NHj,  etc.)  that  physically  cause  the  carbonaceous  melt  to  foam.  If  a  phos¬ 
phoric  acid  is  used,  amides/amines  act  as  dehydrating  agents  promoting  phos¬ 
phorylation  and  enhancing  the  conversion  of  the  carbon  in  the  polyol  to  char. 

In  addition  to  the  gases  from  amides/amines,  the  conventional  intumes- 
cent  system  includes  a  halogenated  organic  material  or  other  suitable  materi¬ 
al.  These  materials  release  a  large  volume  of  nonflammable  gases,  forming 
bubbles  which  result  in  the  characteristic  puff  or  foam  of  an  intnmescent  sys¬ 
tem.  While  water  vapor  and  carbon  oxides  also  may  be  released  during  dehydra¬ 
tion  and  thermal  decomposition  of  the  polyol,  these  gases  will  function  as 
blowing  agents  only  if  they  are  released  at  the  proper  time  (i.e.,  after  melt, 
but  before  aolidification  occurs).  In  order  to  insure  sufficient  gas  release 
at  the  proper  point  in  the  reaction,  blowing  agents  with  predetermined  decom¬ 
position  temperatures  are  employed,  e.g.,  chlorine  containing  paraffins  (HC1 
release  over  a  temperature  range  of  160°  to  350*0  or  melamine  (release  of 
ammonia,  low  molecular  weight  hydrocarbons,  and  carbon  oxides  between  300*  and 
350*0 . 


As  the  reaotion  nears  completion,  gelation  and  finally  solidification 
occurs.  This  solid  is  in  the  form  of  multicellular  foam  which  contributes  to 
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insolation  of  tlxe  substrate.  Typically,  this  char  has  a  volume  that  is  from  5 
to  200  times  the  original  volume  of  the  coating. 

As  can  be  seen  readily  from  the  above  discussion,  if  any  one  of  the 
reactions  does  not  proceed  at  the  required  time,  intumescence  will  not  occur. 
Thus,  the  temperature  at  which  the  specified  reactions  take  place  is  extremely 
important. 

The  binder  insures  a  continuous  and  uniform  structure  throughout  the 
eoating  system  and  helps  to  insure  a  uniform  foam  structure.  It  essentially 
holds  the  intumescent  system  plus  inert  fillers  together.  It  can  be  a  thermo¬ 
plastic  resin,  or  a  convertible  resin  which  hardens  by  either  oxidation  or 
polymerization.  The  resin  film  must  not  be  so  inflexible  that  it  interferes 
with  intumescence.  The  binder  must  also  incorporate  the  expected  properties 
of  a  coating  system: 

(1)  Stable  for  long  times  in  its  pre-application  state  (e.g.,  storage 
in  a  can), 

(2)  Noncorrosive, 

(3)  Easily  applied, 

(4)  Present  an  aesthetically  pleasing  appearance,  and 

(5)  Resistant  to  soiling,  water,  UV,  etc. 

An  intumescent  coating  system  must  generally  be  able  to  provide  all  the  advan¬ 
tages  of  a  conventional  paint  system  plus  the  added  benefits  of  fire  protec¬ 
tion. 
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3.0  ANALYTICAL  MODELING 


3.1  Introduction 


Heat  is  initially  dissipated  by  absorption  into  the  coating  mass  which 
has  a  given  specific  heat  or  heat  capacity.  As  the  heat  capacity  is  exceeded 
during  the  thermal  flux,  the  components  of  the  filler  are  then  activated  by 
various  mechanisms.  A  hydrated  solid  will  first  begin  to  melt  or  liquify, 
dissipating  some  thermal  energy  in  this  process,  then  it  proceeds  to  distill 
away  increments  of  the  bound  water,  further  dissipating  energy  as  it  distills. 
A  state  of  viscoelasticity  is  achieved  by  the  molten  or  liquified  hydrate 
which  traps  the  evolving  water  vapor  and  other  gases  and  commences  to  expand, 
inflate,  grow,  i.e.,  intumesce.  This  also  dissipates  thermal  energy.  Final¬ 
ly,  the  water  or  volatiles  are  gone  and  the  expanding  matrix  releases  the  last 
traces  of  moisture,  becoming  a  solid  anhydrous  mass  of  increased  volume. 

If  the  above  description  occurs  in  the  proper  sequence,  then  the  char 
acts  as  a  thermal  barrier  because  of  its  low  thermal  conductivity.  Ihe  heat 
transfer  to  the  protected  substrate  is,  thus,  significantly  retarded.  This 
matrix  is  also  a  heat-dissipating  medium,  especially  if  the  matrix  expands  to 
a  char  along  with  the  intumescent  fillers  (most  intumescing  systems  have  this 
characterization).  The  resulting  char  is  graphitic  (polyaromatic)  and  is  re¬ 
sistant  to  burning.  Such  chars  are  like  coke  and  charcoal  to  some  degree  in 
that  they  may  glow  and  slowly  oxidize  in  air,  but  do  not  deflagrate  easily. 

As  the  char  glows,  it  begins  to  radiate  back  what  can  be  a  significant  amount 
of  heat. 


An  appropriate  mathematical  model  must  describe  the  various  processes 
discussed  previously.  The  ultimate  goal  of  any  intumescent  system  is  to  pro¬ 
tect  a  substrate,  such  as  the  metal  skin  of  a  warhead.  Indeed,  the  effective¬ 
ness  of  an  intumescing  system  is  based  generally  on  its  ability  to  protect  the 
substrate  thermally.  Typically,  the  effectiveness  is  determined  from  the  tem¬ 
perature-time  history  of  the  substrate.  Thus,  to  compare  the  mathematical 
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description  with  experimental  data,  a  one-dimensional  heat  transfer  model  has 
been  developed  to  describe  the  heat  transfer  through  an  intumescing  system  to 
a  metal  substrate. 


Figure  3-1  depicts  schematically  the  thermal  model  developed.  Figure 
3-l(a)  depicts  the  coating  material  at  time  t  =  0.  At  a  later  time  t^,  heat 
has  caused  the  intumescent  system  to  react,  expanding  and  forming  a  char.  Fig¬ 
ure  3-l(b).  Eventually,  if  the  heat  is  applied  for  enough  time,  the  entire 
coating  will  have  been  transformed  into  char.  The  metal  substrate  is  assumed 
to  have  an  adiabatic  back,  and  the  applicable  equations  are  applied  to  the  in¬ 
tumescent  material  in  differential  elements  dx  from  x  -  0  to  the  outer  sur¬ 
face,  x  =  x_(t).  The  continuity  and  energy  equations,  listed  below,  are  de- 

D 

rived  in  Appendix  A. 


3.2  Continuity  Equations 


dt 


+  P  7T 


3v 

dx 


e 

r 

s 


(3-1) 


dt 


+  P 


3y 
g  dx 


+  f 

8 


(3-2) 


The  density,  p,  represents  the  "solid"  density.  As  the  material  expands 
through  intumescence,  the  total  solid  mass  becomes  less  because  of  outgassing 


of  volatile  products,  f^.  The  mass  lost  to  outgassing  can  be  computed  from 


laboratory  experiments  and  is  a  function  of  the  local  temperature  and  the 


chemical  components  of  the  intumescent  coating.  Ihe  gas  density,  p  ,  is  de- 

8 

pendent  on  the  outgassing  products.  The  outgassing  activity  takes  place  only 


in  the  pyrolysis  zone.  Figure  3-1,  which  is  expanding  at  a  velocity,  v,  due  to 
the  pressure  of  the  outgassed  components  on  the  viscous  "semi-solid"  bulk. 

The  gradient  of  the  velocity  is  non— zero  only  in  the  pyrolysis  zone.  It  is 


zero  in  the  char  region  and  in  the  coating  region  where  f  =  0.  This  is 


equivalent  to  asserting  that  there  is  no  expansion  (or  compression)  of  materi¬ 
al  outside  the  pyrolysis  region. 
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Die  dp/dt  in  Equations  (3-1)  and  (3-2)  refer  to  the  material  deriva¬ 
tive: 


d£ 

dt 


dj> 

at 


+  v 


dx 


(3-3) 


The  reason  for  this  notation  will  he  explained  in  the  discussion  of  the  compu¬ 
tational  algorithms  formulated  to  obtain  numerically  the  time  history  (posi¬ 
tion,  density,  and  temperature)  of  the  intumescent  coating  and  substrate. 


3.3  Energy  Equation 


+  P. 


3T 

dt 


8_  3T 
dx  dx 


v 


3T  _ 

dx  p,lchem 


(3-4) 


In  this  equation,  T  is  the  temperature,  fc  is  the  thermal  conductivity,  C  and 

P 

(C^)g  are  the  solid  and  gas  specific  heats,  h  and  h^  are  the  specific  enthal¬ 
pies  of  the  solid  and  gas,  and  4^^  I*  the  rate  of  heat  per  unit  mass  gener¬ 
ated  by  chemical  reactions.  The  convention  has  been  adopted  where  a  negative 
40hea  tapli**  an  exothermic  chemical  reaction.  The  term  on  the  left-hand  side 
represents  the  sensible  energy  accumulated,  while  the  terms  on  the  right-hand 
side  of  Equation  (3-4)  represent,  respectively,  the  net  heat  conduction  in  the 
solid,  the  heat  absorbed  by  phase  changes,  the  heat  absorbed  due  to  the  expan¬ 
sion  or  convective  process,  and  the  heat  generated  by  chemical  reactions. 
Notice  that  the  various  mechanisms  described  in  subsection  3.1  of  this  section 
are  described  mathematically  in  Equation  (3-4).  Heat  is  conducted  into  the 
material  which  increases  in  temperature.  As  the  temperature  increases,  the 
material  melts,  liquifies,  and  outgases,  described  by  the  second  term  on  the 
right-hand  side,  which  absorbs  energy.  The  expansion  of  the  material,  the 
third  term,  also  absorbs  energy.  The  thermal  conductivity,  k,  can  change 
between  the  virgin  and  char  material,  generally  decreasing  in  value,  the  char 
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being  a  better  insulator.  Additionally,  because  of  the  expansion,  tbe  heat 
must  be  conducted  through  more  material,  thus  further  retarding  the  heat  flow 
to  the  substrate. 


It  has  been  assumed  in  the  derivation  of  these  equations  that  the  gas 
temperature  and  solid,  or  bulk  temperature,  are  identical  at  the  same  posi¬ 
tion,  x.  This  is  not  an  unreasonable  assumption  since  the  phase  transitions 
occur  at  constant  temperature.  Also,  if  the  gas  generated  in  the  interior 
moves  "slow  enough"  to  the  surface,  it  vill  have  essentially  the  same  tempera¬ 
ture  as  the  bulk  material.  Additionally,  the  gas  velocity  has  been  assumed  to 
be  equivalent  to  the  bulk  velocity,  again  a  reasonable  assumption  since  it  is 
the  gas  that  pushes  against  the  "soft  and  pliable"  viscous  material  causing  it 
to  intumesce.  Equation  (3-4)  can  be  slightly  simplified  by  neglecting  the 
pg  (C^)g  terms  since  the  solid  density  is  approximately  a  thousand  times  larg¬ 
er  than  that  of  the  evolved  gases,  while  the  solid  and  gaseous  specific  heats 
are  the  same  order  of  magnitude. 


3.4  Expansion  Factor 


One  additional  equation  is  required  to  relate  the  velocity  to  other 
parameters.  Normally,  a  momentum  balance,  similar  to  the  mass  and  energy  con¬ 
servation  equations,  would  be  performed  to  obtain  the  relationship  between  the 
pressure  generated  by  the  volatile  products,  the  bulk  viscosity  of  the  semi¬ 
fluid  solid,  and  the  velocity.  But  the  pressure  field  and  viscosity,  which 
would  be  strong  functions  of  temperature,  are  not  easily  attainable  experimen¬ 
tally. 


An  alternate  approach  has  been  adopted  analogous  to  an  approach  taken 
by  Clark,  et  al.  [4].  A  characteristic  of  an  intumescent  system  is  that  it 
expands  when  subjected  to  sufficient  heat.  Expansion  factors  of  5  to  10  are 
coanaon,  but  paints  have  been  formulated  with  factors  as  high  as  50  to  60  times 
their  original  thickness.  It  is  not  unreasonable  to  associate  the  amount  of 
expansion  with  the  mass  loss.  That  is,  it  is  the  outgassing  of  pyrolyzed 
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material  that  generates  the  pressure  which  drives  the  expansion.  Bence,  an 
expansion  factor,  E^,  is  defined  to  be  proportional  to  the  quantity  of  out- 

gassed  material : 


(3-5) 


where  at  is  the  original  mass,  m  is  the  final  or  char  mass,  and  m  is  the  cur- 
o  c 

rent  mass.  The  proportionality  is  chosen  such  that  Expression  (3-5)  has  the 
correct  limits: 


The  maximum  expansion  factor  is  given  by  (E.)  .  For  example,  it  is  the  5  to 

X  BtX 

10  just  mentioned. 


The  movement  of  a  "particle"  can  be  written  in  terms  of  its  original 
coordinate  position: 


7 


/ 


Ef Im(x)] 


dz 


(3-7) 


where  y  is  the  coordinate  the  expanding  (or  intumescing)  particle  has  moved 
to,  its  original  location  being  the  coordinate  x.  It  can  be  seen,  by  examin¬ 
ing  Equation  (3-0,  that  if  the  paint  at  a  particular  location  has  not  begun 

to  outgas  (i.e.,  intnmesce),  then  m  =  m  and  the  local  expansion  factor  is  1.0 

o 

[Ej(no>  -  1.0].  Equation  (3-7)  then  implies,  as  it  should,  that  the  particle 

has  not  moved  to  a  new  location,  i.e.,  y  ■  x.  When  the  paint  has  intumesced 

fully,  the  density  is  equal  to  the  char  density,  m  »  m  ,  and  E,(m  )  «  (E,) 

c  to  x  max 

Examination  of  Equation  (3-7)  shows  that  the  position  of  a  particle  is  not 
only  dependent  upon  what  is  occurring  locally,  but  also  depends  upon  the  his¬ 
tory  of  interior  points  of  the  paint.  As  a  "point"  intumesces,  the  expansion 
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factor  associated  with  that  "point"  changes  continuously  from  1.0  to  (E.) 

i  max 

The  ultimate  position  of  a  particle  A  is  the  integral  effect  of  all  the  parti¬ 
cles  interior  to  particle  A. 

The  exponent,  n,  in  Equation  (3-6)  simply  describes  the  dependence  of 
Ef(m)  on  the  change  in  mass  from  mo>  If  the  expansion  occurs  very  early  in 
the  outgassing  stage,  that  is,  small  changes  from  the  original  mass  result  in 
a  significant  portion  of  the  expansion,  then  the  exponent  lies  between  0  and 
1.  If  the  expansion  occurs  only  after  significant  outgassing,  i.e.,  large 
changes  from  the  initial  paint  mass  occur  before  significant  expansion  ensues, 
then  the  exponent  n  is  greater  than  1.  An  exponent  of  1  implies  a  linear 
relationship  between  the  outgassing  products  and  the  expansion.  Figure  3-2 
plots  the  scaled  expansion  factor  versus  the  scaled  mass  for  various  values  of 
the  exponent  n.  The  features  just  described  can  be  observed  by  examining  this 
graph.  For  example,  for  an  exponent  of  10,  when  the  local  density  of  the 
paint  has  undergone  75  percent  of  its  change  from  m^  to  mc,  the  material  has 
expanded  24  percent  of  its  maximnm  value. 

The  particle  velocity,  v,  can  be  found  by  differentiating  Equation 
(3-7)  with  respect  to  time: 


v(x, t) 


iz  = 

dt 


/ 


d(Ef [m(z)]) 

dm 


(3-8) 


After  performing  the  derivative  with  respect  to  m  in  Equation  (3-8)  and  rear¬ 
ranging  terms.  Equation  (3-8)  becomes: 


x 

v(x't)  "  f  [Ef(n)  "  X]  dt  dz  (3"9) 

o 
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Equation  (3-9)  shows  that  the  velocity  at  a  point,  x,  depends  upon  the  accumu- 
lative  effect  of  what  is  occurring  interior  to  that  point.  Even  if  the  local 
change  in  the  coating  mass  (due  to  outgassing)  is  zero,  the  velocity  can  have 
a  non-zero  value  because  of  changes  in  the  mass  at  interior  points. 

3.5  Boundary  Conditions 

The  boundary  conditions  for  the  energy  equation.  Equation  (3-4),  are 
straightforward.  The  condition  at  the  outer  boundary  is  the  requirement  that 
the  heat  flux  be  continuous  at  x  =  x^,  from  which  the  surface  temperature  can 
be  determined: 

.  *  *  (tf  -  Ts)  ♦  •.  •  (•,  V  -  V)  <3-10) 

x "  *B 

where  h  is  the  convective  heat  transfer  coefficient,  e«j  is  the  surface  emis- 
sivity,  Sp  is  the  flame  emissivity,  o  is  the  Stefan-Bol tzmann  constant,  Tp  is 
the  effective  flame  temperature,  and  Tg  is  the  surface  temperature  of  the 
coating,  i.e.,: 

Ts-T(Vt)  (3-11) 

The  symbol  x  **  x^  indicates  that  the  temperature  gradient  is  computed  into  the 
coating  material. 

The  boundary  condition  at  the  interface  of  the  coating  and  substrate, 
x  »  0,  also  requires  that  the  heat  flux  be  continuous  at  the  interface.  How¬ 
ever,  because  of  the  large  value  of  the  thermal  conductivity  for  the  substrate 
(a  steel  plate)  and  the  relative  "thinness"  of  the  plate  (1.6  mm),  essentially 
no  temperature  gradient  exists  through  the  thickness  of  the  plate.  Assuming 
an  adiabatic  back  to  the  substrate  (i.e.,  no  heat  losses  from  the  back  of  the 
plate),  all  heat  crossing  the  boundary  at  x  ■  0  goes  into  heating  the  sub¬ 


strate: 
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P  M  T  ^  = 

V  {  v)p  P  it  dx 

The  subscript  p  refers  to  substrate  (plate)  properties,  with  p,  C  ,  x,  and  T 

P 

referring  to  the  density,  specific  heat,  thickness,  and  temperature  of  the 
plate,  respectively.  This  condition  determines  the  plate  temperature,  vhich 
is  also  the  temperature  at  the  interface  of  the  coating  and  plate: 

T  -  T(0, t)  (3-13) 

P 

The  boundary  condition  on  Equation  (3-1)  is  determined  as  follows.  At 
time  zero,  the  density  at  x  *  is  the  initial  density,  pq: 

P(V°)  ■  "o  (3-14) 

The  surface  temperature,  for  very  early  times,  can  be  estimated  by  solving  the 
heat  equation  for  transient  heat  flow  in  a  semi-infinite  solid  exposed  to  a 
constant  heat  flux,  qQ«  The  surface  temperature,  as  a  function  of  time,  t, 
can  be  computed  from: 


(3-12) 


surface 


T  +  2.  ,1-Y-S-t/jI 
o  k 


(3-15) 


where  Tq  is  the  initial  surface  temperature,  and  a  and  k  are  the  thermal  dif- 

fusivity  and  thermal  conductivity  of  the  coating  material.  For  an  incident 

heat  flux  typical  of  aviation  fuel  fires,  2.7  cal/cm2-s  (10  BTU/ft2-s),  the 

time  for  the  surface  to  reach  800°C  (using  typical  values  for  a  and  k,  and  an 

initial  surface  temperature  of  25'C)  is  computed  from  Equation  (3-15)  to  be  of 

the  order  of  one  second.  Most  intumescent  systems  have  finished  all  activity 

by  800°C.  It  has  just  been  shown  that  the  surface  temperature  surpasses  very 

rapidly  the  temperature  required  to  outgas.  Thus,  at  time  0*,  the  surface 

density  becomes  the  char  density,  p  : 

c 


NADC-822 11-60 


(3-1 


What  mast  be  prescribed  now  is  tbe  form  of  f  ,  which  is  a  function  of 

temperature.  This  couples  Equations  (3-1)  and  (3-4).  f  is  obtained  from 

experimental  data  and  this  procedure  is  discussed  in  Section  3.0.  With  r  a 

$ 

prescribed  function  of  temperature,  we  then  have  the  unknowns,  p,  T,  and  v, 
and  Equations  (3-1),  (3-4),  and  (3-9)  with  appropriate  boundary  and  initial 
conditions. 
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4,0  INTOMES  CENT  PAINT  FORMULATION 


4*1  Formulation 


The  intumescent  paint  formulation  examined  for  this  study  vas  developed 
by  NASA  Ames  Research  Center  (Formulation  EX- 1C- 82 ) ,  and  is  one  that  the  U.  S. 
Navy  has  found  to  be  effective.  It  consists  of  two  parts,  A  and  B,  which  are 
combined  just  before  application.  The  formulation  consists  of  three  major 
components:  the  intumescent  filler  —  borax;  the  binder  —  an  epoxy;  and  a 
fiber  —  glass.  Since  epoxies  tend  to  be  very  hard  and  brittle,  polysulfide 
has  been  added  to  make  the  binder  more  flexible.  Hydrated  salts  such  as  sodi¬ 
um  tetraborate  decahydrate,  or  borax,  function  by  a  different  reaction  mecha¬ 
nism  than  the  one  described  in  Section  2.0.  First,  they  decompose  endothermi¬ 
cally.  This  process  absorbs  a  large  quantity  of  heat  and  effectively  prevents 
the  transfer  of  energy  through  the  coating.  The  second  process  involves  tran- 
spirational  cooling  by  the  water  vapor  produced  from  the  decomposition  as  it 
leaves  the  coating  surface.  The  sublimation  of  water  molecules  cools  the 
upper  layers  of  the  coating  by  convection  and  reduces  the  thermal  transport  of 
energy  through  the  boundary  layer  region  immediately  above  the  surface.  The 
third  process  involves  the  formation  of  a  viscous  material.  Because  of  the 
blowing  action  of  the  water,  the  molten  material  is  transformed  into  an 
expanded  char  which  solidifies  and  provides  thermal  insulation  [5]. 

The  formulations  for  Part  A  and  Part  B  are  given  in  Table  4-1,  along 
with  their  respective  concentration  by  mass  of  the  various  constituents  for 
Part  A,  Part  B,  and  Part  A+B.  Table  4-2  lists  the  percentage  by  weight  of 
each  of  the  constituents.  It  is  these  numbers  which  will  be  used  later  in  the 
modeling  to  combine  the  constituents. 
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Table 


Component 

Poly sulfide 


DKP-30 


EPON  828 


Borax 


Glass 


*  A  and  B  sized  at  a 


4-1.  Intone scent  Paint  Formulation  (Masses) 
NASA  Formulation  EX-1C-82 


ensity  (g/cm  ) 


1.270 

0.973 

1.170 

1.730 

2.500 

TOTAL 


Concentrations  [By  Mass  (g)] 


A  B  (A+B) » 

45.10  -  45.10 

6.50  -  6.50 

-  46.00  46.00 

63.20  45.00  108.20 

0.32  0.25  0.57 

115.12  91.25  206.37 


ratio  (by  mass)  of  1.26  to  1. 
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Component 


Polysulf ide 


DMP-30 


EPON  828 


Borax 


Glass 


Part  A 


Part  B 


*  A  and  B  mixed 


Table  4-2.  Intuaescent  Paint  Formulation 
(Percentage  Concentration) 

NASA  Formulation  EX-1C-82 


Concentrations  (%  By  Yeight) 


B 


39.18 


5.64 


50.41 


54.90 


49.32 


0.28 


0.27 


at  a  ratio  (by  mass)  of  1.26  to  1. 
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4.2  Density  Determination 


4.2.1  Initial.  Balk  Density 


The  balk  density  of  the  coating  system  can  be  computed  from  the  densi¬ 
ties  of  the  constituents.  The  total  mass  of  the  coating  system  can  be 
obtained  by  summing  the  mass  contributions  of  each  of  the  constituents: 

**total  ^borax  +  ^DMP-SO  +  ^polysulfide  +  **EPON  +  Mglass  ^  ^ 
The  bulk  density,  p,  is  defined  by: 


P  V 


(4-2) 


where  V  is  the  total  volume  of  the  coating  system  and  X  is  the  percentage  of 

th  ^ 

total  mass  of  the  j  constituent.  For  example: 


Sor..  "tot.l  *  >W.«  (4-3> 

Assuming  the  total  volume  is  the  sum  of  the  volumes  of  the  constituents  (gen¬ 
erally  a  valid  assumption  except  for  solutions): 


x.  m.  .  , 

j  total 


pj  vj 


»JXJV 


(4-4) 


and  hence: 


P 


J 

£  * 

j-i 


i  pi 


Equation  (4-5)  can  then  be  used  to  compute  the  bulk  density. 


(4-5) 
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As  already  stated,  the  intumescent  paint  formulation  consists  of  two 
parts,  A  and  B,  which  are  combined  jnst  before  application.  Using  the  densi¬ 
ties  and  percentage  by  weight  from  Tables  4-1  and  4-2,  respectively,  the  theo¬ 
retical  balk  densities  of  Part  A,  Part  B,  and  Part  A+B  have  been  computed. 
Table  4-3. 


Table  4-3.  Theoretical  Density  Determinations 

/  3 

g/cm 


Part 

A 

1.509 

Part 

B 

1.450 

Part 

A+B 

1.483 

A  bulk  density  for  Part  A+B  was  determined  by  two  methods.  The  first 
method  involved  the  preparation  of  two  sample  plates.  The  plates  were  weighed 
before  and  after  coating,  the  thickness  of  the  plate  and  plate  plus  coating 
measured,  and  the  cross-sectional  areas  of  the  coatings  determined.  The  den- 

3 

sity  determined  by  this  manner  was  approximately  1.0  g/cm  ,  substantially  less 

3 

than  the  theoretical  value  of  1.483  g/cm  .  Examination  of  the  process  of  mix¬ 
ing  Parts  A  and  B  and  application  to  the  plates  indicates  that  air  gets 
trapped  in  the  volume.  After  the  coating  has  dried,  "air  bubbles"  can  be 
identified.  The  more  viscous  the  formulation  during  mixing  and  application, 
the  more  air  that  can  get  trapped.  The  trapped  air  manifests  itself  in  an 
increased  thickness  of  the  coating,  thereby  reducing  the  experimentally  deter¬ 
mined  bulk  density.  There  is  no  effective  control  (except  experience)  of  how 
much  air  is  entrapped  (with  subsequent  increased  thickness).  When  the  densi- 

3 

ties  for  the  first  two  plates  came  out  to  approximately  1.0  g/cm  ,  a  third 
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plate  vai  prepared.  Great  care  was  taken  to  avoid  air  entrapment,  bat  obser¬ 
vations  daring  application  indicated  the  "air  pockets"  were  still  unavoidable. 

3 

A  value  of  1.22  g/cm  was  determined  for  this  third  plate. 

The  second  method  was  a  specific  gravity  measurement  of  an  oncured  Part 

2 

A+B  sample,  which  yielded  a  density  of  1.49  g/cm  ,  in  very  good  agreement  with 
the  theoretical  limit.  Balk  densities  for  Part  A  and  Part  B  were  also  deter¬ 
mined  through  specific  gravity  measurements  and  foand  to  be  1.32  and  1.31 
g/cm  ,  respectively.  Again,  entrapped  air  increased  the  volume  of  these  sam¬ 
ples. 


The  mass  of  the  applied  coating  is  important  since  enthalpy,  specific 
heat,  etc.,  are  given  in  units  of  energy  per  gram  of  material.  An  effective 
thickness,  which  would  be  the  thickness  if  no  air  were  entrapped  in  the 
applied  coating,  will  be  used  for  modeling  purposes. 

4.2.2  Final,  Char  Density 


The  final  density  is  related  to  the  total  mass  lost  and  the  expansion 
of  the  coating.  The  initial  mass,  mQ,  and  the  final  mass,  m^,  are  given  by: 


o 


p  dx  »  p 
o  *0 


X 

o 


p  dx 
c 


'•Xf 


(4-6) 


(4-7) 


where  Xq  is  the  original  thickness  and  X^  is  the  final,  or  expanded,  thickness 
of  the  coating: 
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X 


f 


Conservation  of  mass  requires; 


m 

o 


+  m 


lost 


(4-8) 


(4-9) 


or,  rearranging  terms: 


But  Equation  (4-9),  making  use  of  Equations  (4-6)  and  (4-7),  can  also  be  writ¬ 
ten  as: 


P  X 
*o  o 


xt  *  *i„,t 


(4-11) 


Thus : 


'•  x»  '  (EfU  x»  +  ■»  1  ' 


“_f 

a 

o 


(4-12) 


Equation  (4-12)  can  be  solved  for  the  char  density: 


(Ef) 


(4-13) 


max 


As  will  be  seen  in  the  next  section,  the  quantity  (mf/mo)  is  determined  from 
thermo gravimetric  analyses.  The  amount  of  expansion,  (EJ  ,  is  determined 

I  BIX 

experimentally  from  measurements  of  the  thickness  of  the  coating  material 
before  and  after  intumescence. 
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5.0  TGA  DATA  —  MASS  LOSS 


5.1  Introduction 


Thermogravimetric  analysis  (TGA)  measures  and  records  the  mass  of  a 
sample  versus  temperature.  A  known  mass,  typically  a  few  milligrams,  is 
placed  on  an  accurate  balance.  The  sample  is  heated  at  a  prescribed  rate 
(e.g.,  20°C/min)  and  the  sample  mass  is  recorded  as  a  function  of  temperature. 

Generally,  the  rate  of  mass  loss  as  a  function  of  temperature  is  also  recorded 
at  the  same  time. 

i 

The  TGA  data  for  the  various  constituents  need  to  be  put  into  a  form 
which  can  be  easily  handled  by  the  numerical  analysis  scheme.  Two  procedures 
will  be  presented,  fitting  to  Arrehnius'  equation  and  Fourier  series  represen¬ 
tation.  The  Fourier  series  representation  was  finally  chosen  as  the  best  pro¬ 
cedure  for  reasons  to  be  discussed  below. 

5.2  Curve  Fitting 

5.2.1  Arrehnius*  Equation  to  Represent  Rate  of  Mass  Loss 

It  is  well  known  that  the  rate  of  many  chemical  reactions  is  strongly 
dependent  upon  the  temperature,  and  Arrehnius  found  that  this  rate,  r,  could 
be  expressed  as: 


r  -  B.-E/RT 


(5-1) 


where  B  is  the  frequency  factor,  E  is  the  activation  energy,  R  is  the  gas  con¬ 
stant,  and  T  is  the  absolute  temperature. 

The  problem  with  Equation  (5-1)  is  that  it  does  not  go  to  zero  when 
the  reactants  have  been  completely  used.  In  the  particular  case  of  a  reactant 
being  one  of  the  constituents  of  an  intumescent  coating  formulation,  then 
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eventually  all  the  "volatiles"  will  be  outgassed.  Equation  (5-1)  can  be  modi¬ 
fied  to  make  it  go  to  zero  by  multiplying  by  a  term  that  looks  like: 


(5-2) 


where  m  is  the  current  mass,  m^  is  the  initial  mass,  and  m^  is  the  residual 
mass.  The  exponent  b  allows  for  nonlinearities. 


The  rate  of  mass  loss  then  can  be  expressed  as: 


(5-3) 


where  B,  E/R,  and  b  have  to  be  determined  by  curve  fits  to  the  data.  One 
method  for  fitting  Equation  (5-3)  to  the  data  is  to  first  take  the  natural 
logarithm: 


Next,  consider  the  change  in  the  above  expression  as  sequential  data  points 
are  analyzed: 


Note  that  In  B  is  a  constant  so  that  the  change  in  In  B  is  zero.  By  consider¬ 
ing  adjacent  data  points  from  the  TGA  curve,  the  change  in  the  natural  log  of 

d(a'aj) 

dT 
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the  change  in  the  inverse  temperature,  and  the  change  in  the  natural  log  of 
the  mass  can  all  be  computed.  By  plotting: 


A  In 


£d(m/m.^/dT  J  A  In  ^m  - 


A  (1/T) 


versus 


A  (1/T) 


(5-6) 


then  the  slope  of  the  straight  line  is  the  exponent  b,  and  the  intercept  is 
-E/R.  Once  b  and  -E/R  are  determined.  Equation  (5-4)  can  be  used  to  determine 
the  frequency  factor,  B. 

An  alternate  method  for  determining  the  three  unknown  coefficients  is 
to  perform  a  least  squares  curve  fit  of  Equation  (5-4),  which  has  the  form: 


where : 


bz 


(5-7) 


(5-8) 


a  *  In  B 


(5-9) 


c  =  E/R 


and: 


where  the  constants  a,  c,  and  b  are  to  be  determined. 


(5-10) 


(5-11) 


One  last  comment  before  further  discussion  of  these  two  curve  fit  pro¬ 
cedures.  The  constituents  of  the  paint  formulation  are  compounds  with  various 
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types  of  chemical  bonds.  These  different  bonds  are  broken  at  different  tem¬ 
peratures.  That  is.  more  energy  (as  indicated  by  a  higher  temperature)  is 
required  to  break  some  bonds  relative  to  others.  Thus,  we  can  expect  that  the 
curve  fits  must  be  applied  only  over  certain  regions  of  the  TGA  data,  and  that 
different  activation  energies,  frequency  factors,  etc.,  will  be  obtained  for 
each  region.  These  parameters  vill  be  unique  to  the  particular  reaction  tak¬ 
ing  place. 

The  procedure  indicated  by  Equation  (5-5)  was  used  to  analyze  the  TGA 
data.  Hovever,  it  produced  unsatisfactory  results.  This  can  be  understood  by 
looking  at  the  circled  portion  of  the  TGA  curve  in  Figure  5-1.  Expression 
(5-2)  is  only  important  as  m  -*■  m^.  Otherwise,  the  exponential. in  Equation 
(5-3)  dominates.  The  precision  and  accuracy  at  vhich  data  can  be  interpolated 
is  at  best  0.025  on  the  ordinate,  and  2.5°C  on  the  abscissa,  which  thus  limits 
the  number  of  data  points  for  curve  fitting.  This  precision  is  just  not  suf¬ 
ficient  to  allow  "good"  curve  fits. 


The  procedure  indicated  by  Equation  (5-7)  was  partially  successful,  but 
the  constants  determined  by  the  curve  fits  are  extremely  sensitive  to  the  val¬ 
ue  used  for  the  initial  mass,  m^,  for  the  onset  of  a  reaction  and  the  final  or 

residual  mass,  m  ,  at  the  termination  of  the  reaction.  As  will  be  discussed 
r 

below,  the  temperature  range  over  vhich  the  various  reactions  take  place  must 
be  determined  graphically.  In  the  procedure  described  below,  we  used  the 
original  mass,  mo>  as  our  reference  point. 

By  plotting  ln[d(m/mo)/dT]  versus  1/T  (the  natural  logarithm  is  auto¬ 
matically  taken  by  plotting  on  semi-log  paper),  the  temperature  range  over 
which  a  particular  reaction  takes  place  can  be  ascertained.  An  example  is 
shown  in  Figure  5-2  for  polysulfide,  where  two  regions  are  depicted  by  the 
symbols  A  and  D.  In  fact,  least-squares  curve  fits  to  the  linear  portions 
give  the  slope,  E/R,  and  the  intercept,  B,  to  the  expression: 


-  Be-E/RT 


(5-12) 
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where  mQ  is  the  original  mass  of  the  intnmescent  system.  The  only  caution 
with  using  Equation  (5-12)  is  that  d(m/mo)/dT  must  be  set  to  zero  when  the 
mass  reaches  the  residual  mass  for  that  particular  reaction. 

Summarizing,  the  Arrehnius  equation  is  quite  adequate  for  describing 
the  mass  loss  of  those  components  which  undergo  a  significant  mass  loss  at 
primarily  one  temperature,  e.g.,  polysulfide,  DMP-30,  and  the  EPON  resin. 
Polysulfide  could  be  represented  by  Equation  (5-12)  by  dividing  the  data  into 
five  temperature  regions.  However,  this  procedure  becomes  more  difficult  when 
a  TGA  curve  is  much  more  complicated.  The  composite  formulation  (A+B) 
required  15  temperature  intervals  [and  15  equations  of  the  form  of  Equation 
(5-12)]  to  describe  the  TGA  data.  The  data  have  to  be  plotted. to  determine 
all  the  different  temperature  regions  for  the  purpose  of  curve  fitting.  In 
the  case  of  Part  A+B,  to  obtain  a  good  curve  fit,  the  physical  significance 
of  activation  energy  has  been  lost  (some  regions  have  a  positive  exponent). 
Thus,  this  procedure  is  not  too  attractive  if  a  number  of  TGA  data  curves  are 
to  be  analyzed. 

5.2.2  Fourier  Series  Representation 

a.  Methodology 


Because  of  the  aforementioned  difficulties  of  using  the  Arrehnius 
equation,  we  have  turned  to  a  Fourier  series  representation  to  describe  the 
mass  loss  as  a  function  of  temperature.  As  will  be  shown,  a  Fourier  series 
representation  has  some  very  nice  features  for  handling  and  combining  the 
data. 


The  mass  loss  as  a  function  of  temperature  can  be  written  as  a  sum 
of  cosine  terms  times  their  respective  Fourier  coefficients: 
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where  T  is  the  temperature  interval  of  interest.  The  Fourier  coefficients 

max 

are  given  by: 


n  =  0,  1,  2,  3  . . . 


(5-14) 


The  Fourier  coefficients  can  be  found  from  the  experimental  data  by  numerical¬ 
ly  performing  the  integration  using  the  trapezoidal  rule: 


a 

n 


f 


fn  n 


cos 


max 


(5-15) 


where  T  =*  T  The  f(T. )'s  are  the  mass  loss  as  a  function  of  tempera- 

JBIAZ  BAX  J 

ture  and  can  be  read  right  from  the  TOA  data.  It  was  found  that  keeping  30 
terms  of  the  Fourier  series  (31  coefficients  counting  aQ)  described  the  exper¬ 
imental  TGA  curve  quite  nicely.  Tables  5-1  and  5-2  list  the  Fourier  coeffi¬ 
cients  for  the  active  components  of  the  intumescent  formulation. 


The  *q  term  represents  the  average  mass  loss  over  the  entire  tem¬ 
perature  range.  The  total  mass  loss  between  any  two  temperatures,  T^  and  T^, 
is  found  by  integrating  Equation  (5-13): 
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Table  5-1.  Coefficients  for  Fourier  Series  Representation 
of  Mass  Loss  Versus  Temperature 


a  x  10  [Equation  (5-13)  With  T  =  830«C] 


max 


Polvsulf ide 


DMP-30 


EPON  Resin 


0 

-23.837200 

-23.851600 

-23.8121000 

1 

-10.201400 

-15.803700 

-  8.0017800 

2 

13.292900 

-  1.518170 

16.1376000 

3 

20.331300 

10.953300 

18.1010000 

4 

5.625640 

18.957100 

-  0.0527777 

5 

-14.034000 

17.616600 

-15.7047000 

6 

-16.288300 

5.643540 

-12.5818000 

7 

1.117050 

-  6.697220 

4.9775600 

8 

15.545200 

-14.096400 

14.7672000 

9 

10.359100 

-14.225300 

5.2917300 

10 

-  6.147620 

-  6.852190 

-  8.3865000 

11 

-14.613400 

3.278050 

-  9.4520600 

12 

-  5.645370 

9.687940 

-  0.0238021 

13 

9.514870 

11.073600 

7.7369200 

14 

12.167500 

6.136110 

5.6203100 

15 

0.457644 

-  1.730880 

-  3.0595700 

16 

-10.492600 

-  7.176350 

-  6.8591900 

17 

-  8.483380 

-  7.788330 

-  1.6365800 

18 

3.139390 

-  4.523490 

4.6420500 

19 

10.039000 

1.554480 

4.4529700 

20 

4.351090 

5.750160 

-  0.9230400 

21 

-  5.796920 

5.758830 

-  4.6691000 

22 

-  8.005040 

2.933670 

-  2.2591100 

23 

-  0.605511 

-  1.176090 

2.8329100 

24 

6.916620 

-  4.526830 

3.7055700 

25 

5.661040 

-  4.444060 

-  0.2895870 

26 

-  2.151600 

-  1.856220 

-  3.3918300 

27 

-  6.608690 

1.365300 

-  1.7943500 

28 

-  2.850380 

3.396200 

1.9160100 

29 

3.727200 

3.188520 

2.7002200 

30 

5.069970 

1.078080 

-  0.1179260 

Borax 

7.947930 

6.200150 

3.056320 

0.770384 

0.421483 

1.076450 

1.812130 

2.229640 

2.172330 

1.829370 

1.235300 

0.659015 

0.391571 

0.678590 

1.238240 

1.705800 

1.956020 

1.914740 

1.667770 

1.243320 

0.668549 

0.156630 

1.164170 

1.997790 

2.408130 

2.319620 

1.902730 

1.356190 

0.886634 

0.533482 

0.149958 
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Table  5-2.  Coefficients  for  Fonrier  Series  Representation 
of  Mass  Loss  Versus  Temperature 


a  x  10^ 
n 


n 

Part  A 

0 

-15.084400 

1 

-  7.643500 

2 

3.218150 

3 

8.231690 

4 

3.923760 

5 

-  2.115430 

6 

-  2.478680 

7 

1.303550 

8 

5.974840 

9 

3.998670 

10 

-  1.438150 

11 

-  6.151490 

12 

-  3.785520 

13 

1.970420 

14 

4.887970 

15 

2.820540 

16 

-  1.930820 

17 

-  3.135440 

18 

-  1.071900 

19 

2.455260 

20 

1.968220 

21 

-  1.253590 

22 

-  3.860280 

23 

-  2.742830 

24 

0.529961 

25 

2.202990 

26 

1.571850 

27 

-  0.643321 

28 

-  0.91  83  50 

29 

0.429189 

30 

2.225970 

[Equation  (5-13)  With  T  = 

max 


Part  B 

-15.2637000 

-  5.7946500 
7.3794600 
7.8223500 
0.7606680 

-  4.1354100 

-  2.9925600 
2.9898400 
5.9568400 
1.9698300 

-  3.3808300 

-  4.6454300 

-  1.1002200 
2.9765600 
1.9932400 

-  1.5676600 

-  1.6635700 
0.6137730 
1.5838800 
0.5327720 

-  1.6746400 

-  2.4188700 

-  0.5342380 
1.1583100 
0.9595690 

0.1457600 

0.0311966 

1.0346300 

1.8587800 

0.8723340 

-  0.7435300 


830#C] 


Part  A  +  B 

-15.8770000 

-  3.3232700 
2.6975300 
3.1181500 
1.7166100 

-  1.4819600 

-  2.0689700 
2.8357900 
2.9845200 
1.3247800 

0.3722910 

-  2.7066400 

-  0.7185560 
1.0548600 
1.4344500 

0.5010040 

0.2879580 

-  0.3001990 
1.0024100 
1.6266700 

1.1428100 

-  0.6697590 

-  1.1008800 

-  0.5433120 

-  0.5706260 

0.0991600 

-  0.2825600 

-  0.9209250 

-  0.1361500 
0.3600440 

0.1014060 


38- 


NADC-822 11-60 


(5-17) 


The  total  mass  loss  from  T  =  0  to  a  temperature  T  is  obtained  directly  from 
Equation  (5-16)  by  letting  =  0  and  -*■  T: 


(v). 


lost 
0  ->  T 


a  T 


o 


2 


(5-18) 


or  the  mass  remaining  at  the  temperature  T  is: 


-fr' 


lost 
0  -*•  T 


(5-19) 


Figures  5-3  and  5-4  depict  m/mo  versus  the  temperature  for  the 
Fourier  series  representation  and  the  experimental  data  for  each  of  the  con¬ 
stituents  of  the  formulation.  As  can  be  seen,  the  Fourier  series  representa¬ 
tions  provide  excellent  reproductions  of  the  TGA  data.  Figures  5-5  and  5-6 
depict  d(m/mo)/dT  versus  temperature  [from  Equation  (5-13)]  for  each  of  the 
constituents  of  the  formulation.  As  will  be  explained  in  Section  5.3,  these 

data  are  used  for  computing  f  [Equations  (3-1)  and  (3-4)].  It  is  interesting 

8 

to  note,  from  Figure  5-5,  that  even  though  the  magnitude  of  d(m/mQ)/dT  for 
borax  is  substantially  smaller  than  the  other  three  components,  borax  is  con¬ 
sidered  to  be  the  "blowing  agent."  The  borax  curve  is  repeated,  with  a  dif¬ 
ferent  ordinate  scale  in  Figure  5-6.  It  is  also  worthwhile  to  note  that  a 
substantial  amount  of  activity  is  still  occurring  in  Part  A+B,  the  actual  in- 
tumescent  formulation,  for  temperatures  in  excess  of  400  to  500°C.  This  ac¬ 
tivity  is  not  present  in  any  of  the  other  figures  and  is  probably  the  conse¬ 
quence  of  the  chemical  bonding  between  the  two  elements  of  the  epoxy  binder. 
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(a)  Poly sulfide  (b)  DMP-30 


(c)  EPON  828  Resiu  (d)  Borax 


Figure  5-3.  Fourier  Series  Representation  of  TGA  Data  for 
Polysulfide,  DMP-30,  EPON  Resin,  and  Borax 
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TEf*>  <C) 

(a)  Polysulfide 


• ea  in  2aa  aa  sat  aae  sea  7»a  aa*  » a  laaa 


TEMP  <C> 

(b)  DMP-30 


(c)  EPON  828  Resin 


(d)  Borax 


Figure  5-5.  Fourier  Series  Representation  of  d(m/mo)/dT  for 
Polysulfide,  DMP-30,  EPON  Resin,  and  Borax 


-42- 
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The  different  constituents  can  be  mixed  by  simply  summing  the  Four¬ 
ier  coefficients  times  their  respective  relative  initial  concentrations  by 
weight,  Xj : 


d/m/m 


«  J  ,  ^ 

-“composite  “"r  J 


E  *»  5j  o#,(t; 


(5-20) 


ft) 


lost 

0  -*■  T 


•  S  ft). 


composite 


'lost 
0  T 


The  mass  versus  temperature  curve  can  now  be  computed  for  the  mixture  using 

Equations  (5-19)  and  (5-21).  This  procedure  has  been  followed  to  form  the 

(m/m  )  versus  temperature  curves  for  Part  A,  Part  B,  and  Part  A+B.  Figure  5-7 
o 

compares  the  simple  "mixing  theory"  for  Part  A  and  Part  B.  The  relative  con¬ 
centrations  of  the  individual  constituents  are  given  in  Table  4-2.  Overall, 
the  agreement  is  quite  good.  The  d(m/mo)/dT  versus  temperature  curves  corre¬ 
sponding  to  the  curves  in  Figure  5-7  are  shown  in  Figure  5-8. 


c-*»ei*  la/wi 
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TEMP  <C> 

(a)  Part  A 


TEMP  <C) 

(b)  Part  B 


Figure  5-8.  Summation  of  the  Various  Constituents,  by  Concentration 
to  Form  Part  A  and  Part  B  d(m/mo)/dT  Versus  Temperature  Curves 
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Likewise,  the  entire  formal ation.  Part  A+B,  can  be  approximated  as 
a  snm  of  its  parts.  Figure  5-9.  The  agreement  between  the  simple  mixing  the¬ 
ory  and  the  data  is  not  as  good  as  for  Part  A  or  Part  B.  However,  considering 
that  no  chemical  bonding  has  been  accounted  for,  which  particularly  occurs 
when  the  two  components  of  the  epoxy  axe  mixed  together,  the  overall  trend  of 
the  experimental  data  is  given  quite  nicely  by  the  simple  mixing  theory.  Fig¬ 
ure  5-10  depicts  the  curves  for  the  corresponding  rate  of  mass  loss  versus 
temperature.  What  remains  to  be  determined  is  if  the  model's  heat  transfer 
predictions  differ  significantly  for  Part  A+B  compared  to  Part  A+B  from  the 
mixing  theory. 


The  relative  concentration  of  a  particular  constituent  is  tempera¬ 
ture  dependent,  and  thus  must  be  computed  at  each  temperature.  Since  there 
exists  a  temperature  gradient  through  the  intunescent  coating,  the  concentra¬ 
tion  of  the  various  constituents  is  a  function  of  depth  into  the  insulating 

system.  The  local  constituent  concentration,  c.,  is  computed  from: 

J 


V"  *T 


("'•o),  (*«), 


(5-22) 


where  (m/m  ).  is  the  local  residual  mass  remaining  of  constituent  j.  The 

®  J  ^ 

original  mass  of  the  j  constituent  is  given  by  the  total  mass,  m  ,  times  its 

o 

original  concentration,  X. : 

J 


(a  \  =  m  X 

'  %  01 


(5-23) 


Thus,  the  local  constituent  concentration  is  given  by: 


). x; 


-o'j  -j 


(5-24) 
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TEMP  (C) 

(a)  (Part  A)  •  XA  +  (Part  B)  •  XB  -*•  Part  A+B 


TEMP  <C> 

(b)  Summation  of  Constituents  (Borax,  etc.)  by 
Concentration  to  Form  Part  A+B 


Figure  5-9.  Summation  of  the  Components,  by  Concentration,  to 
Form  Part  A+B  Mass  Loss  Versus  Temperature  Curves 
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TEMP  CC> 

a)  (Part  A)  •  Xa  +  (Part  B)  •  Xfi  Part  A+B 


TEMP  CO 

(b)  Summation  of  Constituents  (Borax,  etc.)  by 
Concentration  to  Form  Part  A+B 


Figure  5-10.  Summation  of  the  Components,  by  Concentration,  to 
Form  Part  A+B  d(m/mo)/dT  Versus  Temperature  Curves 
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It  becomes  evident  that  c.  is  a  function  of  the  temperature  since  (m/m  ).  is 

J  o  j 

temperature  dependent.  Note  that  (m/mQ)j  is  straightforward  to  compute  from 
Equation  (5-19). 


5.3  Interpretation  of  Data  for  Analytic  Model 


The  quantity  f  in  Equation  (3-1)  has  units  of  mass  per  unit  volume  per 

unit  time.  We  must  perform  the  appropriate  conversions  on  the  TGA  data  to 

transform  the  data  into  a  form  compatible  with  the  definition  of  f  .  The 

g 

change  in  mass  with  respect  to  time  can  be  written  as: 


dm  =  dm  dT  _  d(m/mo)  dT 
dt  "  dT  dt  “o  dT  dt 


(5-25) 


The  sweep  rate  of  the  temperature,  dT/dt,  is  selected  by  a  switch  on  the  in¬ 
strument.  For  the  TGA  runs  we  have  analyzed,  a  sweep  rate  of  20°C  per  minute 
was  used,  thus  dT/dt  =  1/3  #C/s,  but  in  general,  for  small  increments  of  tem¬ 
perature  and  time: 


dT  AT 
dt  “  At 


(5-26) 


The  change  in  density  of  the  coating  is  due  to  mass  loss  and  the  expan¬ 
sion  of  the  coating.  The  expansion  of  the  coating  is  essentially  one-dimen¬ 
sional,  so  a  differential  volume  element,  AV,  is  related  to  its  initial  volume 
(AV)o  by: 


AV  -  A  (Ax)  «AE.  (Ax)  =  E.  (AV)  (5-27) 

X  o  I  o 

where  A  is  the  cross-sectional  area,  E.  is  the  expansion  factor,  and  (Ax)  is 

x  o 

the  original  thickness  of  the  element. 


The  mass  lost  by  outgassing  is  given  by: 


-50- 


X 
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6.0  DSC  DATA 


6.1  Introduction 


A  differential  scanning  calorimeter  (DSC)  measures  and  records  the 
energy  necessary  to  maintain  thermal  equilibrium  between  a  test  sample  and  a 
chosen  reference.  The  sample  of  known  mass  (generally  a  few  milligrams)  is 
heated  at  a  constant  temperature  rate  (typically  20°C/min)  and  the  rate  at 
which  energy  is  absorbed  (endothermic)  or  emitted  (exothermic)  is  recorded 
versus  temperature.  Specific  heats  are  thus  proportional  to  the  deviation  of 
the  sample  thermogram  from  the  baseline.  A  phase  change  (solid  to  liquid, 
liquid  to  vapor,  or  solid  to  vapor)  is  represented  by  a  peak.  The  area  of  the 
peak  is  proportional  to  the  heat  of  fusion,  vaporization,  or  sublimation.  The 
DSC  curve  for  the  intumescent  formulation  (Part  A+B)  is  shown  in  Figure  6-1. 

To  quantify  the  data  from  the  instrument,  three  standards  were  run 

(indium,  tin,  and  zinc)  which  have  sharp  endothermic  reactions  of  156. 6°C, 

231. 9°C,  and  419. 4°C,  respectively,  with  enthalpy  changes  of  0.9657,  1.154, 

and  1.7226  cal/g.  Since  the  masses,  m  (mg),  of  each  of  the  samples  were 

2 

known,  and  the  areas,  A  (in.  ),  of  the-  endothermic  peaks  were  measured  care¬ 
fully  by  a  planimeter,  three  calibration  constants,  e  ,  were  obtained  using: 

c 


e 

c 


AH  •  m 

A  •  B  •  Aq 
^r 


(6-1) 


where  AH  is  the  enthalpy  change  (mcal/g),  B  is  the  time  base  setting  (sec/ 

in.),  and  Aq^  is  the  y-axis  range  setting  of  the  instrument  (mcal/sec/in. ) . 

The  three  calibration  constants  were  then  combined  into  an  overall  calibration 

constant,  E,  as  a  function  of  temperature  by  fitting  the  e  ' s  to  a  parabolic 

c 

curve: 


E 


2.024  x  10 


-6 


T2  + 


1.714  x  10"3  T  +  0.6476 


(6-2) 


where  T  is  in  degrees  Celsius 
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We  note  here  that  it  is  important  that  the  sample  should  be  tested  in 
an  inert  atmosphere,  such  as  argon,  to  avoid  exothermic  reactions  resulting 
from  the  oxidation  of  long-chain  polymers. 


Enthalpy  changes.  Ah  (cal/g),  were  computed  from: 


Ah  =  “  (B  •  E  •  Aq^  (6-3) 

2 

where  A  is  the  area  of  the  peak  (in.  )  [measured  carefully  by  using  a  planime- 
ter],  B  is  the  time  base  setting  (sec/in.),  E  is  from  Equation  (6-2),  and  Aq^ 
is  the  y— axis  range  (mcal/sec-in. ) .  The  change  in  sample  mass.  Am  (mg),  is 
computed  from  the  TGA  data.  [Note  that  it  has  been  assumed,  and  not  unreason 
ably,  that  the  enthalpy  change  is  associated  with  a  phase  change  of  solid  or 
liquid  to  a  gas  so  that  there  is  mass  loss.  This  assumption  is  consistent 
with  the  7GA  data.  Solid-to-liquid  phase  changes  could  be  accounted  for,  but 
examination  of  TGA  data  generally  shows  mass  loss  during  a  phase  transition. 
This  assumption  will  be  discussed  further  in  the  section  on  reproduction  of 
the  DSC  data.]  The  enthalpies  are  tabulated  in  Table  6-1. 


6.3  Specific  Heats 


Specific  heats,  (cal/g-°C), 


were  determined  from: 


(6-4) 


where  H  is  the  heating  rate  (deg/sec),  AT  is  the  difference  in  y-axis  deflec¬ 
tion  between  the  sample  and  reference  background  curves  at  the  temperature  of 
interest  (in.),  m  (mg)  is  the  same  mass  [which  is  a  function  of  temperature 
and  is  obtained  from  the  TGA  data],  Aq^  is  the  y-axis  range  (mcal/sec-in.), 
and  E  is  from  Equation  (6-2).  Tables  6-2  and  6-3  tabulate  the  specific  heats 


Table  6-1.  Enthalpy  Changes 
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Table  6-2.  Specific  Heats  of  Individual  Constituents 


Temperature 


Spe c ific  Heat  ( cal/g-* 


Polysulfide  DMP-30  EPON  Resin 


0.347 

0.372 

0.386 

0.454 

0.602 

0.376 


0.362 

0.487 

0.180 

0.180 


0.180 

0.180 

0.180 

0.180 


0.180 

0.180 

0.180 

0.180 


0.180 

0.180 

0.180 

0.180 


0.180 


Borax 


Glass  Fiber 
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Table  6-3.  Specific  Heat  of  Mixed  Components 


Temperature  (°C) 


mil 


Part  A 


Specific  Heat  (cal/gm-°C 
Part  B 


5 


T 


A 
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as  a  function  of  temperature.  If  an  exothermic  reaction  takes  place,  then  the 
procedure  adopted  was  to  keep  a  constant  specific  heat  and  balance  the  energy 
by  a  chemical  reaction  term,  discussed  in  the  next  subsection. 

6.4  Heat  Generation  Term 


After  analyzing  the.  DSC  data,  we  concluded  that  some  exothermic  reac¬ 
tions  were  taking  place  at  elevated  temperatures.  To  properly  account  for 
exothermic  chemical  reactions,  an  extra  term  had  to  be  included  in  the  energy 
equation.  This  term  was  (-p  4^^)  •  where  4^^,  *■*  th®  rate  of  heat  generated 
per  unit  mass  by  chemical  reactions  (e.g.,  cal/gm— s).  The  convention  has  been 
used  were  a  negative  4^^  implies  heat  generation. 


For  those  materials  where  the  DSC  curve  goes  above  the  background,  an 


energy  balance  was  used  to  compute  4 


chem* 


4  v  +  C 

chem  p 


AT  =  /E  '  Aqr\  AY  AT 
At  **  l  H  J  m  At 


(6-5) 


where  AT/At  is  the  heating  rate  (e.g.,  1/3  °C/s),  H  is  the  heating  rate  (deg/ 

sec)  [note  that  H  =  AT/At),  E  is  the  cell  calibration  coefficient,  Aq^  is  the 

y-axis  range,  AY  is  the  difference  in  y-axis  deflection  between  the  sample  and 

reference  background,  and  m  is  the  residual  sample  mass  (obtained  from  the  TGA 

curve).  The  specific  heat,  C  ,  was  assumed  to  be  constant  for  all  tempera- 

P 

tures  once  the  DSC  curve  crossed  above  the  background.  The  4  .  *s,  as  a 

cnem 

function  of  temperature,  are  given  as  Tables  6-4  and  6-5  for  each  of  the  com¬ 
ponents  and  the  mixtures. 


6.5  Comments  on  the  Accuracy  of  Tabular  Values 


Specific  heats  and  4cjlea  values  for  some  of  the  constituents  (polysul- 
fide,  DMP-30,  and  EPON  Resin  828)  are  quite  suspect  at  elevated  temperatures. 
For  temperatures  hiaher  than  300  to  400°C.  the  remainina  mass  of  material  is 
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less  than  a  few  percent  of  the  initial  mass.  With  snch  a  small  residual  mass, 

the  computed  C  's  and  q  v  '*  are  extremely  sensitive  to  the  stability  and 
p  cneffl 

accuracy  of  the  DSC  apparatus.  The  measured  differences  between  the  sample 
run  and  the  reference  background  are  probably  in  the  "noise"  range  of  the  par¬ 
ticular  instrument  used.  There  exists  a  much  higher  confidence  in  the  com¬ 
puted  C  'a  and  &  .  's  for  those  substances  where  the  residual  mass  is  at 
p  ^chem 

least  30  percent  of  the  initial  mass. 


To  solve  numerically  the  differential  equations  which  apply  at  all 
locations  within  the  intumescent  coating,  the  coating  material  is  divided  into 
increments.  Figure  7~1.  Because  the  model  is  one-dimensional,  each  increment 
can  be  represented  by  a  grid  point.  The  differential  equations  are  then  writ¬ 
ten  as  difference  equations  and  applied  at  each  of  the  grid  points.  Thus, 
each  of  the  variables,  e.g.,  T,  p,  and  v,  is  computed  at  each  of  the  grid 
points.  The  scheme  we  have  chosen  is  to  let  the  variable  be  representative  of 
the  region  (zone)  to  the  right  of  the  grid  line.  For  example, .the  values  for 
temperature  and  density  for  the  cross-hatched  zone  in  Figure  7-1  are  given  by 
T^  and  where  the  subscript  k  refers  to  the  ktlx  grid  point.  Thus,  T^  and  p^ 
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are  the  temperature  and  density  for  the  region  between  k*Ax  <.  x  <  (k+l)-Ax, 

where  k  runs  from  zero  to  kmax.  With  this  convention,  it  is  easily  seen  that 

the  values  for  T.  and  p,  represent  boundary  conditions  since  no  intumes- 
rmax  xmax 

cent  material  lies  to  the  right  of  the  kmax  grid  point. 

The  overall  approach  is  to  move  forward  in  time  by  computing  all  the 
variables  at  a  new  time  step  in  terms  of  known  values  from  the  previous  time 
step,  and  values  already  computed  at  the  current  time  step.  This  explicit 
formulation,  as  it  is  called,  is  straightforward,  but  the  time  step  must  be 
within  certain  bounds  for  the  numerical  procedure  to  remain  stable. 

A  unique  feature  is  the  manner  in  which  the  equations  are  solved.  The 
energy  equation  is  solved  in  an  Eulerian  (fixed)  coordinate  system  to  deter¬ 
mine  the  temperature  field  in  order  to  compute  the  mass  loss  (which  is  temper¬ 
ature  dependent).  After  computing  the  mass  loss,  the  energy  equation  is  cor¬ 
rected  (by  subtracting  the  convective  term  and  updating  the  pyrolysis  term), 
and  then  the  continuity  and  energy  equations  are  solved  in  a  Lagrangian  (mate¬ 
rial)  coordinate  system.  The  Lagrangian  coordinate  system  is  the  "natural" 
coordinate  system  because  the  7GA  data  are  Lagrangian  in  nature  (since  the 
intumescent  system  expands  during  outgassing).  However,  a  Lagrangian  coordi¬ 
nate  system  was  not  adopted  for  the  overall  computational  approach  because  the 
accuracy  of  the  computational  scheme  would  decrease  with  increasing  spatial 
increments.  A  velocity  is  associated  with  each  grid  point,  and  the  grid 

expands  according  to  the  velocity  field  and  the  time  increment  (e.g..  Ax  = 

move 

v.At).  The  expanded  grid  is  then  rezoned  back  to  its  original  mesh.  Because 
the  coating  system  expands,  the  total  thickness  of  the  computational  mesh 
grows  larger.  New  zones  are  added  when  appropriate. 

7.2  Logic  Flow 

The  logic  flow  is  presented  now  to  give  the  reader  an  overall  view  of 
the  structure  of  the  numerical  algorithms.  Descriptions  of  the  numerical  rep¬ 
resentations  of  the  equations  are  given  in  detail  in  Appendix  B. 


NADC-822 11-60 


Hie  procedure  to  compute  the  temperature,  velocity,  and  density  at  the 
grid  points  is: 

(1)  Compute  the  temperature  at  all  grid  points  by  integrating  Equation 
(3-4). 

(2)  From  the  temperatures,  compute  the  mass  loss  in  each  zone,  i.e., 
outgassing.  Equation  (5-30). 

(3)  From  the  mass  loss,  compute  the  velocity  at  all  grid  points.  Equa¬ 
tion  (3-8) . 

(4)  Compute  the  density  at  all  grid  points  by  integrating  Equation 
(3-1)  —  the  decrease  in  bulk  density  at  each  point  is  the  result 
of  losing  mass  and  expansion. 

(5)  Compute  the  expansion  of  the  grid. 

(6)  Rezone  the  grid  points  back  to  their  original  locations  —  add  grid 
points  as  necessary  to  account  for  expansion. 

This  process  is  repeated  for  each  time  step  until  program  termination.  Logic 
flow  diagrams  depicting  the  logic  flow  of  the  numerical  procedure  are  given  in 
Appendix  B. 
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8.0  MODEL  CALCULATIONS 

Numerical  solutions  have  been  obtained  for  the  analytic  model  presented 
in  the  previous  sections.  The  results  of  the  different  problems  or  parametric 
studies  vill  be  presented  in  the  following  subsections. 

8.1  Reproduction  of  DSC  Data 


Vith  the  information  from  the  tables  entered  into  the  computer  program, 
a  computer  run  was  made  to  simulate  a  DSC  test.  Many  things  must  work  cor¬ 
rectly,  in  concert,  for  the  computational  model  to  reproduce  the  actual  DSC 
experimental  curve,  including  the  mass  loss,  search  subroutines,  enthalpy 
changes,  specific  heat,  chemical  energy  generation,  and  expansion  of  the  intu- 
mescent  coating.  The  total  remaining  mass  is  important  because  of  the  speci¬ 
fic  heat  and  chemical  energy  generation  terms.  Likewise,  the  rate-of-mass 
loss  is  important  since  the  total  heat  absorbed  during  a  phase  transition  is 
the  enthalpy  change  times  the  mass  loss.  The  expansion  of  the  coating  must 
occur  correctly  or  else  the  computed  densities  will  be  in  error,  and  when  in¬ 
tegrated  over  the  thickness  of  the  expanded  coating,  will  result  in  an  incor¬ 
rect  total  remaining  mass. 

Figure  8-1  reproduces  the  data  from  the  DSC  experiment  for  Part  A+B. 
Also  plotted  are  the  predictions  from  the  analytic  model.  The  curves  agree 
extremely  well  in  regions  where  there  are  no  phase  changes.  The  two  enthalpy 
peaks  are  present  in  the  computational  predictions,  but  do  not  exactly  re¬ 
produce  the  data.  The  model  assumption  that  the  total  enthalpy  change  was  due 
to  a  phase  change  involving  mass  loss,  must  be  violated  to  some  extent.  It 
would  be  extremely  difficult  to  differentiate  the  quantity  of  heat  going  into 
solid-to-liquid  phase  changes  versus  solid-to-gas  or  liquid-to-gas  (with  a 
corresponding  mass  loss).  However,  the  areas  of  the  two  peaks  have  been  re¬ 
produced,  which  implies  that  the  total  heat  balance  has  been  faithfully  repli¬ 
cated.  Thus,  the  fact  that  the  enthalpy  peaks  have  not  been  exactly  repro¬ 
duced  is  not  of  great  concern.  Therefore,  since  the  current  model  correctly 
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Figure  8-1.  Model  Comparison  With  Experimental  DSC  Data 
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predicts  the  overall  features  of  the  phase  transition,  it  is  felt  that  the 
present  modeling  assumptions  are  adequate.  Solid-to— 1 iquid  phase  changes 
could  be  incorporated  into  the  model  at  a  later  date  if  desired. 

8.2  Composite  DSC  Data 


Section  5.2.2(b)  showed  the  results  of  mass  loss  versus  temperature 
using  the  assumption  that  the  sum  of  the  constituents  might  approximately 
equal  the  whole,  i.e.,  the  actual  formulation  of  Part  A+B  (Figure  5-9).  The 
same  can  be  done  from  an  energy  absorption/generation  standpoint. 

The  DSC  data  in  Figure  8-1  can  be  displayed  in  an  alternate  form  by 
dividing  the  energy  absorbed/generated  by  the  instantaneous  mass  per  unit 
length  (from  the  TGA  data).  This  rate  at  which  energy  is  absorbed  or  gener¬ 
ated  per  gram  of  coating  material  (per  square  centimeter  of  surface  area)  can 
then  be  displayed  versus  temperature.  The  solid  curve  in  Figure  8-2  is  the 
curve  for  Part  A+B.  A  hypothetical  DSC  experiment  was  then  performed  with 
the  computer  model  using  the  concentrations  from  Table  4-2  to  mix  Part  A  with 
Part  B,  and  mix  the  constituents  (polysnlf ide,  borax,  etc.)  to  form  Part  A+B. 
These  results  are  also  shown  in  Figure  8-2.  The  densities  of  the  composite 
formulations  were  determined  from  Equation  (4-5).  It  is  important,  for  model¬ 
ing  considerations,  that  the  total  mass  of  the  system  be  represented  accurate¬ 
ly.  Thermodynamic  quantities  such  as  specific  heat  are  given  in  units,  for 
example,  of  calories  per  degree  Celcius  per  gram.  If  an  alternate  procedure 
is  used  to  compute  the  composite  density,  then  the  coating  thickness  may  have 
to  be  adjnsted  to  insure  that  the  coating  mass  is  modeled  correctly. 

The  mixing  of  the  constituents  follow  the  Part  A+B  curve  much  better 
than  the  mixing  of  Part  A  with  Part  B.  Part  A  plus  Part  B  has  substantially 
larger  endothermic  reactions  at  approximately  75°C  and  150°C.  Evidently,  the 
mixing  of  components  into  Part  A  and  into  Part  B  is  slightly  exothermic  — 
with  the  application  of  heat,  the  breaking  of  chemical  bonds  absorbs  signifi¬ 
cantly  more  energy. 
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Figure  8-2.  Comparison  of  Rate  of  Heat  Absorption  Versus  Temperature 
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Part  A+B  become  exothermic  at  temperatures  in  excess  of  385°C.  How¬ 
ever,  the  unit  used  to  perform  the  DSC  experiment  was  limited  to  temperatures 
less  than  approximately  600°C.  All  heat  capacities  were  extrapolated  for  tem¬ 
peratures  in  excess  of  600°C  by  assuming  the  specific  heat  a  constant  for 
higher  temperatures.  If  a  sample  was  exothermic  at  600°C,  then  the  exothermic 
reaction  was  set  to  zero  at  750°C.  From  a  computational  viewpoint,  this  im¬ 
plies  that  the  exothermic  heat  generation  goes  to  zero  linearly  with  tempera¬ 
tures  from  6 00°C  to  750#C  [see  Equation  (B-52)].  The  "jitter"  in  the  curve 
for  the  constituents  in  the  temperature  range  475°C  to  650°C  can  be  attributed 

to  the  accuracy  of  some  of  the  specific  heat  and  4  v  values.  Tables  6-2  and 

CH61H 

6-4.  Because  the  mass  of  a  sample  becomes  very  small  (e.g.,  polysulfide)  for 
temperatures  in  excess  of  450°C,  small  errors  in  subtracting  the  DSC  data  from 
the  background  can  result  in  significant  changes  in  calculated  specific  heats 
and  4chem's.  The  jitter  is  probably  an  indication  of  the  stability  (accuracy) 
of  the  DSC  apparatus  in  running  a  sample  and  background  at  different  times. 

It  is  evident  from  Figure  8-2  that  the  actual  -ixing  of  the  components, 
with  the  resulting  formulation  of  chemical  bonds,  has  a  significant  influence 
on  the  rate  of  energy  absorption.  We  can  expect  different  temperature-time 
histories  of  the  substrate  when  the  model  is  exercised  using  the  thermodynamic 
data  for  Part  A+B  versus  the  thermodynamic  data  for  the  components  that  make 
up  Part  A+B. 

8.3  Coated  Substrate 


A  steel  sample  plate  (approximately  7.6  cm  x  7.6  cm  x  0.0813  cm)  was 
coated  with  0.0965  cm  of  intumescent  formulation.  It  was  then  inserted  in  the 
furnace  test  facility  at  the  Naval  Air  Development  Center.  The  temperature  of 
the  substrate,  determined  by  a  thermocouple,  was  monitored  versus  time.  The 
expanded  coating  was  measured  after  the  test  and  an  expansion  factor  of  4.47 
was  determined. 
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The  computer  model  was  set  up  with  five  zones  initially,  each  zone 

being  0.0193  cm  thick  for  a  total  thickness  of  0.0965  cm.  A  flame  temperature 

-4 

of  1227°K  (1750°F)  and  a  convective  heat  transfer  coefficient  of  2.81  z  10 
2 

cal/cm  -s-°C  were  used.  The  flame  emissivity  was  assumed  to  be  1.0,  and  the 

surface  emissivity  0.8.  Vith  these  parameters,  the  initial  incident  heat  flux 
2  2 

was  2.68  cal/cm  -s  (9.89  BTU/ft  -s)  with  radiation  providing  90  percent  of  the 
heat  flux.  This  heat  flux  approximately  reproduces  the  conditions  of  the  test 
chamber. 


The  first  parameter  investigated  was  the  effect  of  the  exponent  in  the 
expansion  factor  [Equation  (3-6)].  For  values  of  n  less  than  1.0,  the  expan¬ 
sion  occurs  early  in  the  outgassing  stage.  For  values  of  n  greater  than  1.0, 
the  expansion  occurs  only  after  significant  outgassing  (Figure  3-2).  As  can 
be  seen  in  Figure  8-3,  the  time  at  which  the  expansion  begins  has  a  rather 
dramatic  effect  on  the  temperature  of  the  substrate. 

The  next  parameter  investigated  was  the  exponent,  q,  in  the  assumed 

transition  law  for  the  thermal  conductivity  [Equation  (B-54)].  (The  values 

-3 

chosen  for  thermal  conductivities  were:  virgin  state,  k  «  J.5  x  10  cal / 

-3 

cm-s-#C,  and  char  state,  k^  3  2,0  x  10  cal/cm-s-°C. )  The  larger  the  expo¬ 
nent,  the  faster  k  approaches  k  for  a  given  density,  the  logic  being  that 
o  c 

the  char  may  be  essentially  formed  during  the  early  stages  of  mass  loss.  Fig¬ 
ure  8-4  presents  the  effects  of  the  exponent  in  the  thermal  conductivity  law 
for  the  case  of  n  =  0.6  from  Figure  8-3.  Except  for  early  times  where  the 
substrate  heats  more  slowly  for  larger  values  of  q,  the  curves  are  essentially 
identical.  Thus,  the  effect  of  the  exponent  q  is  to  delay  heating  of  the  sub¬ 
strate  for  early  times  by  retarding  the  heat  flow  with  a  smaller  thermal  con¬ 
ductivity. 

Figure  8-5  presents  the  influence  of  the  maximum  expansion  of  the  coat¬ 
ing  system  on  the  thermal  response  of  the  plate.  With  the  assumption  that  the 
char  thermal  conductivity  is  relatively  independent  of  density  (the  density 
decreases  with  expansion),  the  thickness  of  the  char  layer  has  a  dramatic 
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influence  on  the  substrate's  thermal  response.  It  is  interesting  to  note  here 
that  the  measured  expansion  was  found  to  be  4.47  times  its  original  thickness. 
A  thick  sample  of  the  coating  (1.0  cm)  was  placed  on  a  substrate  and  exposed 
to  a  heat  source  until  it  had  completely  intumesced.  This  coating  expanded  to 
2.6  cm.  Hie  next  section  will  briefly  discuss  this  observation. 

Plotted  in  each  of  the  figures  are  the  experimental  data  for  the  plate 
test.  The  model  predictions  do  not  exactly  replicate  the  experimental  data. 
Model  assumptions  and  deficiencies  are  discussed  in  the  next  section,  but  some 
observations  and  conclusions  can  be  drawn  from  the  model's  predictions  never¬ 
theless. 


The  experimental  data  begins  rising  rather  quickly,  but  the  slope 
decreases  for  100°C  <  T  <  155°C,  after  which  the  rate  of  heating  increases  but 
is  less  than  the  initial  heating  rate.  The  model  predicts  this  "bend"  in  the 
data,  after  which  the  heating  rate  becomes  steeper.  This  bend  in  the  data  and 
model  predictions  occurs  when  the  temperature  of  the  coating  in  contact  with 
the  substrate  reaches  the  second  endothermic  peak  shown  in  the  DSC  data  (see 
Figure  8-1).  Note  that  if  the  heat  flux  to  the  substrate  is  sufficiently 
large.  Figures  8-3  and  8-3,  the  energy  absorbed  by  the  endothermic  reaction  is 
relatively  insignificant  with  respect  to  the  total  heat  reaching  the  substrate 
and  the  bend  does  not  appear.  The  first  endothermic  peak,  because  it  is 
smaller  and  because  the  thermal  conductivity  of  the  virgin  material  is  larger 
than  that  of  the  char,  is  more  difficult  to  see  in  the  model's  predictions  but 
it  is  discernible.  The  resolution  of  the  experimental  data  is  not  adequate  to 
discern  this  first  endothermic  peak.  However,  possibly  for  reasons  discussed 
in  the  next  section,  the  model  does  not  retard  the  incoming  heat  flux  suffi¬ 
ciently  to  reproduce  the  extent  of  the  bending  evident  in  the  experimental 
data. 


From  Figure  8-3,  it  is  evident  that  the  model  does  indicate,  for  better 
agreement  with  experimental  data,  that  the  expansion  should  occur  early  in  the 
mass  loss  process  (n  <  1).  This  is  in  excellent  agreement,  from  a  ohemistry 
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point  of  view,  that  the  blowing  agent  zesnlts  from  the  hydrating  of  water  from 
borax.  Examination  of  TGA  and  DSC  data  shows  that  this  occurs  in  the  early 
stages  of  the  mass  loss. 

Figure  8-5  might,  at  first  glance,  seem  to  indicate  that  the  largest 

expansion  factor  [(E.)  =  6.0]  yields  closer  agreement  with  the  experimental 

i  max 

data.  However,  for  n  =  0.6,  it  underpredicts  the  initial  heating  rate  and 
then  has  "too  much"  curvature  (i.e.,  is  insulating  too  well)  at  late  times. 

The  rate  of  heating  of  the  substrate  for  temperatures  greater  than 
175°C  is  approximately  replicated  by  the  model.  Figure  8-6,  for  any  of  the 
curves  where  (Er)  =  4.5  and  n  <  0.  Very  minor  differences  exist  in  the 

curves  in  Figure  8—3  for  n  =  0  or  Figure  8-4  (q  =  1,  2,  or  3)  for  temperatures 
greater  than  175°C.  This  implies  that  the  heat  flow  is  being  controlled  by 
the  charred  region  of  the  coating.  The  curve  in  Figure  8-6  was  obtained  by 
shifting  the  model  predictions  in  time  to  match  the  experimental  data.  This 
would  indicate  that  late  time  heat  flow  is  being  treated  properly.  The  model 
predicts  a  little  more  curvature  than  the  experimental  data,  but  then  the  mod¬ 
el  has  assumed  the  rear  surface  of  the  substrate  is  perfectly  insulated.  Any 
heat  loss  from  the  back  surface  of  the  experimental  plate,  even  though  it  may 
be  well-insulated,  would  tend  to  "straighten"  the  experimental  data  curve. 

In  summary,  the  model  appears  to  predict  very  early  heating  rates  and 
the  very  late  time  heating  rates  reasonably  accurately  when  the  expansion  pro¬ 
cess  occurs  early  in  the  mass  loss  process  (0.6  <  n  <  0.8),  for  a  maximum 
expansion  factor  of  approximately  4.5.  The  model  computes  a  heating  rate  a 
little  too  high  in  the  intermediate  times,  but  does  predict  a  relatively  sharp 
increase  in  the  heating  rate  at  approximately  155°C.  For  best  overall  agree¬ 
ment,  the  model  calculations  indicate  that  (E.)  should  be  approximately 

i  ntz 

4.5,  in  excellent  agreement  with  the  experimentally  determined  expansion  fac¬ 
tor  of  4.47. 
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Figure  8-6.  Comparison  of  Late  Time  Heating  Rate  of  Substrate 

With  Model  Predcitions 


NADC-822 11-60 


9.0  CONCLUSIONS 

A  thermodynamic  heat  transfer  model  has  been  formulated  for  intumescent 
systems.  This  model  includes  mass  loss  and  expansion  of  the  coating  in  com¬ 
puting  the  temperature  history  of  a  metal  substrate.  It  should  be  pointed  out 
that  the  temperature-time  history  of  the  substrate  is  a  function  of  the  heat 
flux  applied  at  the  surface  of  the  coating  and  the  chemistry  and  physics  with¬ 
in  the  intumescing  material.  Experimental  data  exists  only  for  the  substrate. 
To  compare  the  model's  predictions  with  experimental  data  requires  the  appli¬ 
cation*  of  a  heat  flux  at  the  boundary  and  computing  the  heat  transfer  to  the 
substrate.  For  the  computer  calculations  presented  in  the  previous  section, 
approximately  18,000  time  increments  (cycles)  were  required  for  each  computer 
solution.  If  the  physics  is  not  modeled  correctly,  substantial  errors  can  be 
expected  since  an  error  would  be  "propagated"  from  the  surface  to  the  sub¬ 
strate,  accumulating  in  time.  The  effects  of  several  parameters  have  been 
shown.  These  parameters  are  describing,  in  a  sense,  something  about  the  phys¬ 
ics  and  chemistry  of  the  intumescing  and  heat  transfer  processes. 

There  exist  uncertainties  in  the  thermal  conductivities  of  the  virgin 

and  char  states,  though  the  model  computations  would  seem  to  indicate  that  we 

have  reasonable  estimates.  The  actual  transition  of  k  to  k  has  been  assumed 

v  c 

to  be  a  function  of  the  density  change.  Figure  8-4  showed  that  there  is  some 
dependence  in  the  heating  rate  for  different  values  of  the  exponent  in  Equa¬ 
tion  (B-54),  but  alternate  forms  of  "transition"  functions  could  be  explored. 

Thermodynamic  data  are  known  only  to  600°C,  whereas,  the  mass  loss  data 
are  known  as  a  function  of  temperature  to  830°C.  But,  the  temperatures  pre¬ 
sent  in  the  fire  environment  will  cause  the  intumescent  material  to  reach  tem¬ 
peratures  in  excess  of  900°C  to  1000°C.  Data  should  be  obtained  for  these 
higher  temperatures. 

The  major  assumption  in  the  model  is  the  assumed  form  of  the  expansion 
law.  Currently,  the  expansion  law  states  that  the  expansion  depends  only  upon 
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the  change  in  mass  relative  to  the  total  change,  or  in  a  more  general  sense, 
is  dependent  only  on  the  mass  loss: 


E 


f 


More  generally,  the  expansion  factor  probably  is  a  function  of  not  only  the 
mass  loss,  but  the  rate  of  mass  loss,  and  the  temperature: 


Since  the  foaming  and  intumescence  of  the  coating  depends  upon' the  generation 
of  gases  to  push  the  material  out,  it  stands  to  reason  that  the  rate  of  vapor 
or  gas  generation  would  control  the  expansion  process.  Also,  the  material 
must  be  fluid  enough  (or  soft  enough)  to  be  pushed  or  expanded,  which  results 
in  the  temperature  dependence.  The  dependence  of  the  model  on  alternate  forms 
of  the  expansion  factor  should  be  investigated.  Alternate,  and  perhaps  more 
realistic  forms,  for  the  expansion  law  could  account  for  the  discrepancy  in 
the  model's  predictions  and  the  experimental  data.  For  example,  it  has  al¬ 
ready  been  indicated  from  exercising  the  model  that  the  expansion  should  occur 
early  in  the  mass  loss  process.  Perhaps  the  expansion  should  all  occur  over  a 
narrow  temperature  range  (which  would  correspond  to  a  small  portion  of  the 
mass  loss  curve).  Experimental  observations  have  shown  that  very  thick  coat¬ 
ings  do  not  expand  to  the  same  extent  as  thinner  applications.  A  temperature 
and  rate  of  mass  loss  dependence  would  imply  that  the  temperature  gradient 
within  the  intumescing  material  could  be  important  in  the  expansion  process. 


An  analytic  model  for  heat  transfer  through  an  intumescing  system  has 
been  presented.  It  has  been  shown  that  it  can  predict  features  observable  in 
experiments,  and  explain  the  physical/chemical  processes  causing  these  fea¬ 
tures.  The  model  has  some  shortcomings,  but  areas  of  improvement  have  been 
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delineated.  As  the  model  continues  to  develop,  it  will  quantify  our  under¬ 
standing  of  certain  features  of  intumescing  systems,  and  perhaps  become  a  use 
fnl,  predictive  tool  for  formulating  new  or  improved  intumescent  coatings. 
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APPENDIX  A 

DERIVATION  OF  CONTINUITY  AND  ENERGY  EQUATIONS 


A-l 
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A.l  Conti natty  Equations 

The  density  of  the  intnmescent  system  decreases  due  to  expansion  and 
the  loss  of  mass  from  pyrolysis.  Figure  A-l  depicts  a  region  with  thickness 
6x  and  cross-sectional  area  A.  Let  the  mass  in  the  region  be  characterized  by 
a  mass  density,  p.  The  rate  of  change  of  solid  (bulk)  mass  contained  in  the 
region,  d(p  6x  A)/dt,  must  equal  the  net  mass  flow  across  the  boundaries  of 
the  closed  region.  The  mass  flow  into  the  region  is  given  by  p(x,  t) • v(x, t) »A 
where  p(x, t)  is  the  bulk  density  and  v(x,t)  is  its  velocity  at  the  coordi¬ 
nate  position,  x.  The  mass  flow  out  of  the  region  is  given  by  -p(x+6x, t) • 
v(x+4x, t).A  where  the  minus  sign  denotes  mass  leaving  the  region.  Note  that 
the  density  and  velocity  are  computed  at  the  coordinate  position  of  the  second 
boundary,  x+Sx.  If  the  temperature  is  high  enough,  pyrolysis  occurs  and  gas 
is  generated.  This  gas  generation  is  at  the  expense  of  the  mass  in  the  bulk 

material  and  is  designated  by  -  f  [T(x>]  Sx  A.  f  6x  A  is  the  rate,  then,  at 

B  B 

which  gaseous  mass  is  generated,  and  it  is  shown  to  be  an  explicit  function  of 
the  temperature,  T. 


Thus,  conservation  of  mass  implies: 
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~  (p  6x  A)  *  p(x,t)  v(x,t)  A  - 

p(x+6x, t)  v(x+8x, t)  A  -  f  [T(x)l  6x  A  (A-l) 

The  second  term  on  the  right-hand  side  of  Equation  (A-l)  can  be  expanded: 

p(x+8x)  v(x+8x)  A  = 

£p.(x,t)  +  8x  +  ...  J  £v(x,t)  +  8x  +  ...  J  A  « 

p(x,t)  v(x,t)  A  +  p(x,t)  6x  A  + 

v(x,t)  8x  A  +  o[(6x)2]  (A-2) 


Thus.  Equation  (A-l)  becomes,  as  the  infinitesimal  element  8x  -*>  0: 


tTU)1 


(A-3) 


A  similar  equation  can  be  written  for  the  change  in  the  mass  of  the  gas 
due  to  outgassing  during  pyrolysis: 


It  *  ’  sf * ",  “  -  * f, 11(1,1 


(A-4) 


where  the  subscript  g  denotes  gas.  The  action  of  pyrolysis  is  the  source  for 
changing  the  gas  density.  It  was  assumed  in  the  derivation  of  Equation  (A-4) 
that  the  gas  velocity,  v  ,  is  equivalent  to  the  bulh  velocity,  v.  It  is  the 
gas  which  "pushes"  against  the  molten  and/or  pliable  viscous  material  causing 
the  coating  to  intumesce.  Thus,  it  is  reasonable  that  the  two  velocities  are 


A-4 
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approximately  equal.  It  is  true  that  the  gas  eventually  diffuses  to  the  sur¬ 
face  as  the  expanding  coating  begins  to  harden  and  form  the  char  residue,  but 
little  error  should  be  incurred  by  the  approximation  of  the  tvo  velocities  be¬ 
ing  equal. 

A. 2  Energy  Equation 


The  energy  equation  is  derived  similarly  to  the  continuity  equations. 
Figure  A-2  depicts  a  fixed  region.  The  rate  of  change  of  energy  within  the 
region  equals  the  net  energy  flow  into  the  region.  The  change  in  energy  is 
given  by  d(ph  6x  A)/dt,  where  h  is  the  specific  enthalpy.  This  change  in 
energy  results  from  heat  conduction,  and  the  mass  flow  of  enexgy  (both  bulk 
and  gaseous  material)  across  the  boundaries: 
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Figure  A-2.  Energy  Conservation 
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where  the  eroes-sectionel  tree  has  been  divided  out  of  each  term.  Performing 
''he  Taylor  seriea  expansion  about  (x, t),  cancelling  terms,  and  letting  the  in¬ 
finitesimal  thickness  &x  go  to  zero  results  in: 


a 


kill  ♦  iliSl  +  fi  .  2_  k  21  .  i_ 

a*  rH.i..  i,  a.  a. 


a 


at 


che: 


.  ■  al  k  n  -  £  <“rt)  -  al  (V\)  <*'« 


The  terms  on  the  left-hand  side  of  Equation  (A-6)  can  be  expanded: 


ifcA)  +  lifikl  m  h  +  p  l^  +  hl£  +  pli 

at  at  g  at  pg  at  at  p  at 


(A-7) 


Likewise,  the  last  two  terms  on  the  right-hand  side  can  be  expinded: 


iCgT.%1  .  h  LCezl  +  pv  M 


ax 


dx 


ax 


(A-8) 


ah 


dx 


<  9* 


P  ▼  T~^ 
I  3* 


(A-9) 


Inserting  Equations  (A-7)  through  (A-9)  into  Equation  (A-6),  and  utilizing  the 
mesa  conservation  equations.  Equations  (A-3)  and  (A-4),  yields: 


ah 


ah 


",  it  *  *  8t 


•f“k“-(h  -  h)  f  -  vp  -  vp  -r-*  -  p4  - 
dx  dx  \  g  /  g  dx  |  dx  chew 


3h 


(A-10) 


We  now  make  use  of  a  thermodynamic  identity  for  enthalpy: 


dh-C  dT 


u 

+  Vd / 


(A-ll) 
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where  C  is  the  specific  heat  capacity  at  constant  pressure,  T  is  the  tempera- 
P 

ture,  V  is  the  specific  volume,  and  P  is  the  pressure.  The  change  in  pres¬ 
sure,  dP,  is  negligible  in  comparison  to  the  changes  in  temperature  in  the 

intumescent  system.  Also,  the  gas  temperature,  T  ,  can  be  taken  to  be  equal 

8 

approximately  to  the  bulk  temperature  since  (1)  phase  transitions  occur  at 
constant  temperature,  and  (2)  if  the  gas  generated  in  the  interior  moves  "slow 
enough"  to  the  surface,  it  will  have  essentially  the  same  temperature  as  the 
bulk  material.  Thus,  Equation  (A-10)  becomes,  after  insertion  of  Equation 
(A-ll) , : 


chem 


(A-12) 
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B.l  Discretizing  the  Differential  Bqnations 


To  solve  numerically  the  differential  equations  which  apply  at  all 
locations  within  the  intumescent  coating,  the  coating  material  is  divided  into 
increments.  Because  the  model  is  one-dimensional,  each  increment  can  be  rep¬ 
resented  by  a  grid  point.  The  differential  equations  are  then  written  as  dif¬ 
ference  equations  and  applied  at  each  of  the  grid  points.  Thus,  each  of  the 

variables,  e.g.,  T,  p,  and  v,  is  computed  at  each  of  the  grid  points. 

The  overall  approach  is  to  move  forward  in  time  by  computing  all  the 
variables  at  a  new  time  step  in  terms  of  known  values  from  the  previous  time 
step,  and  values  already  computed  at  the  current  time  step.  This  explicit 
formulation,  as  it  is  called,  is  straightforward,  but  the  time  step  must  be 

within  certain  bounds  for  the  numerical  procedure  to  remain  stable. 

B.2  Logic  Flow 


The  computer  program  consists  of  a  main  program,  called  PAINT,  which 
controls  the  flow  of  the  numerical  procedure  by  calling  subroutines.  Each 
subroutine  generally  performs  one  specific  computational  procedure.  Figure 
B~1  presents  the  overall  logic  flow  of  PAINT.  The  subroutines  called  from 
PAINT  may  in  turn  call  a  series  of  subroutines,  each  performing  their  own  spe¬ 
cific  computational  task.  Figure  B-2  presents  subroutine  READ  with  the  list 
of  input  values  required  by  the  program.  This  subroutine  also  performs  vari¬ 
ous  initialization  procedures. 

The  subroutine  which  computes  the  time  step  (or  increment)  is  divided 
into  two  completely  separate  portions  with  their  own  entry  points,  TIMSTP 
(time  step)  and  TIMCHK  (time  check).  Figure  B-3.  A  time  step  is  computed  from 
the  requirements  for  numerical  stability  for  the  energy  equation.  This  opera¬ 
tion  is  performed  in  TIMSTP.  The  continuity  equation  also  has  an  associated 
stability  criterion.  After  computing  the  velocities  and  densities,  the  time 
step  criterion  is  checked  against  the  time  step  currently  being  used.  If  the 
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Figure  B-l.  Flow  Diagram  for  Overall  Logic  Flow 
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Figure  B-3.  Flow  Diagram  for  Time  Step  Determination  and  Stability  Check 
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stability  criterion  is  violated,  the  code  corrects  the  time  step  and  recom¬ 
putes  temperatures,  etc.,  with  the  adjusted  time  step. 

Figure  B-4  presents  the  logic  flow  of  the  energy  equation  for  computing 
the  temperatures  throughout  the  coatputationsl  grid.  Also  displayed  in  Figure 
B-4  is  the  subroutine  CORRECT  which  is  called  after  the  amount  of  outgassing 
is  computed.  This  subroutine  prepares  the  energy  equation  for  Lagrangian  com¬ 
putation  by  subtracting  out  the  convective  term.  Also,  the  subroutine  updates 
the  pyrolysis  term  with  the  information  for  the  current  time  step. 

Figure  B— 5  depicts  the  logic  flow  for  the  continuity  equation,  which 
computes  the  velocities  and  densities  at  the  grid  points.  With  the  velocities 
known,  the  expansion  of  the  grid  is  computed. 

After  expansion  of  the  grid,  it  is  rezoned  back  to  the  original  mesh 
geometry.  The  logic  flow  for  rezoning  is  displayed  in  Figure  B-6.  Upon  com¬ 
pletion  of  rezoning,  all  variables  are  updated  and  the  procedure  is  started 
over  to  compute  the  variables  at  the  next  time  step. 

B.3  Difference  Equations  and  Computational  Procedures 

B.3.1  Introduction 


The  difference  equations  will  be  presented  and  briefly  discussed  since 
some  unique  features  are  present  in  the  numerical  formulation.  The  equations 
are  generally  discussed  in  the  order  in  which  they  are  computed  within  the 

program.  The  following  convention  will  be  used:  a  subscript  k  will  refer  to 

th  i 

the  k  grid  point  and  a  superscript  i  will  refer  to  the  time  step,  e.g.,  T. 

til  i  * 

is  the  temperature  at  the  k  grid  point  at  the  time,  given  by  t  ,  found  by 

summing  the  time  increments  from  t»0: 
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Ill*  derivatives  are  replaced  by  differences,  e.g.,: 


3T  ^  AT 
dx  Ax 


d%  _  AT 
3t  At 


^  "  Xk  “  Xk-1 


At  -  ti+1  -  t1 


(B-2) 


(B-3) 

(B-4) 

(B-5) 


In  general,  all  the  zone  sizes  are  equal  except  for  the  zone  at  the  boundary 
which  varies  in  size.  This  zone,  where  explicitly  necessary,  will  be  desig- 

“t*4  **  %  '  *km.x  '  ‘ta.r-1- 

B.3.2  Energy  Eqnation  and  Its  Boundary  Conditions 


the  energy  equation  consists  of  five  tens.  Each  of  these  tens  will 
be  considered  separately. 

a.  Sensible  Heat 


PC 


8T 

p  at 


(B-6) 


b.  Conduction 


The  difference  fon  for  the  conduction  term  can  best  be  visualized 
by  referencing  Figure  7-1.  First,  the  spatial  derivatives  between  the  k+1  and 
k,  and  k  and  k-1  points  are  considered: 
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then,  the  derivative  between  these  two  points  is  taken: 


3*  -  ax  1/2  RAx),  +  (Ax), 


The  energy  transported  by  conduction,  E1 ,  is  thus  given  by: 


(Ax).  +  (Ax) 


Except  at  the  (kmax-1)  grid  point,  (Ax)t  «  (Ax)2  -  Ax  which  simplifies  Equa 
tion  (B-9).  At  the  (kmax-1)  grid  point,  (Ax)2  ■  Ax^  and  (Ax^  »  Ax. 


c.  firrolysis 


d.  Convection 


32  '  -  f,  '  (f,)‘ 
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e.  Chemical  Eneri 


E.  “  -p6  . 
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f .  Computation  of  Tempera tore 


It  should  ho  noted  that  C  ,  h,  q  .  ,  and  f  all  depend  upon  the 

p  cnem  g 

temperature  and  are  found  in  tables,  or  in  the  case  of  f^,  by  an  analytic  ex¬ 
pression.  The  temperature  at  the  new  time  step  can  no*  be  determined  in  terms 
of  quantities  from  the  previous  time  step: 

C  ■  Tk  ♦  T77V  [E1  'VV  h]  (B-14> 

(cp)t  L  J 


g.  Outer  Boundary  Condition 


The  condition  at  the  boundary  is  given  by  Equation  (3-10)  but  modi¬ 
fied  to  incorporate  heat  absorption/generation  from  phase  changes,  convection, 
and  chemical  reactions.  Equation  (3-10)  is  the  true  boundary  condition  to  the 
differential  equation,  but  ve  shall  account  for  these  other  heat  terms  because 
of  the  finite  size  of  the  grid,  (le  shall  ignore  heat  capacity  effects  in  the 
boundary  condition  because  of  its  detrimental  influence  on  the  computational 
time  step.)  Equation  (3-10)  then  becomes: 


(T.  -  T.-i) 


s-1 


Ax. 


(E2  *  E3  +  AX_ 


Q 


(B-15) 


where : 


T^1 

km  ax 

(B-16) 

T^1 

kmax-1 

(B-17) 

■  S  »  (*»  TJ  -  Tt) 

(B-18) 

Notice  that  in  the  limit  of  Ax^  -*■  0,  the  second  term  in  Equation  (B-15)  goes 
to  zero  and  Equation  (B-15)  reduces  to  Equation  (3-10). 
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Iliis  is  a  nonlinear  equation  in  the  unknown  T# . 

find  T  is  to  compute  a  correction  6T  and  then  use: 
s  s 


The  procedure  to 


jnew  ^  „old 


T  +  6T 

s  s 


(B-19) 


to  compute  a  new  correction.  This  procedure  is  repeated  until  8Ts  reaches 
some  reasonably  small  value: 

°  (T.  '  T.-l)  *  «  *  (E2  *  E3  *  E<)i_1 


+  H  +  4e  8_  <x  T3 
s  F  s 


(B-20) 


h.  Inner  Boundary  Condition 


The  inner  boundary  condition.  Equation  (3-13),  has  a  similar  modi¬ 
fication  to  account  for  energy  absorption/generation  because  of  the  finite 
width  of  the  zone: 


Solving  for  T. 


i+1 


At 

„i+l  -i  . 

>  ( 

ft  -  ft)/4*  +  (EJ  +  E3  *  E4)„  41 

40  *0 

’(PV).t..!  +  "»  Wo  ^  ’ 

(B-22) 


Again  note  that  in  the  limit  of  Ax  -*■  0,  Equation  (B-22)  reduces  to  Equation 
(3-13). 


B-15 
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B.3.3  Mass  Loss 


The  methodology  for  computing  d(m/mo)/dT  is  given  by  Equation  (5-20) 

where  the  temperature  in  Equation  (5-20)  is  given  for  each  zone,  T*+*.  Equa- 

i+1  * 

tion  (5-28)  is  used  then  to  compute  (r  ).  .  The  total  mass  lost  for  each 

I  1 

constituent  is  also  computed  from  Equations  (5-18)  and  (5-19)  so  the  local 
concentrations  of  each  of  the  constituents,  given  by  Equation  (5-22),  can  be 
coaiputed. 


Because  the  Fourier  series  representation  can  "oscillate,"  and  thus 
d(m/mo)/dT  >  0,  a  check  is  put  into  the  algorithm  to  insure  that  mass  is  not 
added  during  a  time  step.  The  computational  technique  is  general  in  that  it 
will  allow  mass  addition  with  a  resultant  decrease  in  total  thickness  (since 
the  resultant  boundary  velocity  will  be  negative).  But  this  is  non-physical. 
Thus,  whenever  Equation  (5-20)  predicts  that  d(m/m  )/dT  is  greater  than  zero, 
d(m/mo)/dT  is  redefined  to  be  zero,  i.e.,  no  suss  loss  is  allowed  during  that 
particular  time  step.  However,  this  procedure  could  ultimately  result  in  too 
much  mass  being  subtracted,  so  another  check  is  put  into  the  code: 


dT 


0  if 


£  a 

i-0  * 


(At)  Ax  < 


ft) 


lost 
0  -  T 


(B-23) 


Jk 


where  the  right-hand  side  of  the  inequality  is  given  by  Equation  (5-21). 
B.3.4  Correct  Energy  Equation 


Equation  (B-14)  is  used  to  compute  the  temperatures  at  the  fixed  grid 
points.  These  temperatures  are  then  used  to  coaqiute  the  amount  of  outgassing. 
However,  the  material  expands  due  to  outgassing.  Thus,  the  Eulerian  computa¬ 
tion  for  temperature  must  be  corrected  back  to  Lagrangian  coordinates.  This 
is  easily  accomplished  since  the  convective  term.  Equation  (B-12)  is  basically 
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the  one  term  that  distinguishes  the  two  coordinate  systems.  In  addition,  the 
pyrolysis  term  is  updated  for  the  enthalpy  change  at  the  current  time  step. 
Hence,  the  temperature  is  redefined: 


(B-24) 


where  E 2  is  given  by  Equation  (B-ll)  and  E^  is  given  by  Equation  (B-12). 

Equation  (B-24)  is  modified  appropriately  [see  Equation  (B-22)]  for  the 
temperature  of  the  substrate,  if1.  The  temperature  at  the  surface  is  recom¬ 
puted  using  the  methodology  of  Section  B. 2. 2(g),  but  modified  to  not  include 
the.  cony.e  c t ive_term  E^.  but  using  jfche.  .updated  pyrolysis  term. 

B.3.5  Expansion  Factor 


The  expansion  factor.  Equation  (3-6),  is  not  difficult  to  calculate. 

The  expansion  factor  essentially  provides  the  Lagrangian  description  of  the 
location  of  a  material  particle  and  that  location  depends  upon  the  past  his¬ 
tory  (i.e.,  the  amount  of  material  outgassed).  The  original  mass  in  a  zone 

(m  ) . ,  the  final  or  char  mass  (m  ).,  and  the  current  zone  mass  m,  are  required 
ok  c  k  th  k 

for  the  computation.  The  mass  lost  in  the  k  zone  is  given  by: 


'  (M> " 

Hence,  the  expansion  factor  is  given  by: 

wr  ■ 1  *  [(m.„  -  *] 

Note  that  the  original  mass  in  a  fixed  zone  changes  with  time  [(mo)^].  Be¬ 
cause  of  expansion,  the  original  mass  will  be  partitioned  into  several  zones 


\i+l 

‘lo,,)k 


(Ml 


-  m 


(B-26) 
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at  a  later  time,  This  tracking  of  masses  is  accomplished  during  the  rezoning 
procedure. 

In  the  local  time  coordinates*  the  system  has  already  undergone  a  cer¬ 
tain  percentage  of  its  total  expansion.  The  local  time  step  expansion  factor 
is  given  by  the  ratio  of  the  expansion  factor  due  to  the  total  mass  loss  at 
the  current  time  step  to  the  expansion  factor  computed  from  the  previous  time 
step: 


(B-27) 


B.3 .6  Continuity  Equation 


a.  Determination  of  Velocity 

The  continuity  equation*  Equation  (3-1)*  is  solved  in  two  steps. 
First,  the  velocity  is  computed  from  the  mass  loss.  The  local  expansion  fac- 
tor.  [<Ef)t  Jloc.j  i.  given  by  Equation  (B-27).  Using  rectangular  quadrature 
for  the  integral  in  Equation  (3-8)*  the  velocity  is  calculated  from: 


b.  Computation  of  Density 


Ax 


(B-28) 


Once  the  velocity  has  been  determined.  Equation  (3-1)  can  be  writ¬ 


ten  as  a  difference  equation: 
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i+1 

k+1 


Ax 


(B-29) 


and  thus,  the  density  at  the  new  time  step  is  given  by: 

>r  -  [^  (or  a*  ]/[ 

Caution  mast  be  exercised,  however,  as  the  amount  of  mass  in  a  zone 
approaches  the  terminal  or  char  mass.  Computationally,  if  the  time  step  is 
too  large,  an  excess  of  mass  conld  be  deleted  from  a  zone.  The  procedure  fol¬ 
lowed  is  discussed  in  the  next  subsection. 

c.  Maas  Iterate 

Because  of  the  finite  time  step,  an  excess  of  mass  could  be  deleted 
from  a  zone  (because  of  outgassing)  which  would  result  in  the  density  falling 
below  the  char  density.  The  algorithm  developed  computes  the  mass  loss  neces¬ 
sary  for  the  zone  to  reach  its  final  or  char  state,  taking  into  account  that 
the  zone  will  expand  during  outgassing  (whenever  the  zonal  mass  is  within  10 
percent  of  the  char  mass):. 

-  4 4*  •  [  (Ef)‘+ii  o. <®-n> 

L  J local 


The  mass  loss  computed  from  the  TEA  data  is  given  by: 


-  (f  V+1 


At  Ax 


(B-32) 


This  quantity  is  compared  against  the  left-hand  side  of  Equation  (B-31).  If 
Expression  (B-32)  is  greater  than  Am^^,  then  (r^)***'  is  reduced  accordingly: 


B-19 
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( Am.  V+1 

/■  \i+l  V  l08Vk  ..  _/a  \H 

Vg)k  (At)  (Ax)  (  g/k 


At  Ax  > 


(B-33) 


This  is  the  rate  of  density  change  due  to  ontgassing  of  volatiles  which  is  re¬ 
defined  by  Equation  (B-33). 


To  compute  (Am^  #)*  .  a  correction  &m  is  computed  from: 


^♦K^r-roci 

.  _  L  J  local 

.  I  . 


p  Ax 
c 


-  1  -  p  Ax 
rc 


ic-r  wi 


(B-34) 


where  [(E^ ^k+^lOCal  giv«a  by  Equation  (B-27).  Then  (A®lofs)k+^  i*  updated 


(4"lo.,)‘+1  *  (i“lo..)‘+1  +  6*t 


(B-35) 


This  procedure  is  iterated  until  bm^  approaches  "zero." 

B.3.7  Movement  of  Grid  Points 

The  velocity  at  each  grid  point  has  been  computed  by  Equation  (B-28) . 
The  material  then  expands  according  to  the  velocity  field,  and  the  expanded 
grid  coordinates  are  given  by: 


i+1  _  i  .  l+l  .  * 
xk  -  *k  +  \  At 


(B-36) 


B.3.8  Rezonini 


At  the  conclusion  of  the  Lagrangian  computations,  the  solution  has  been 
advanced  in  time  with  appropriate  movement  of  the  grid.  The  next  step  in  the 
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procedure  is  to  rezone  or  remesh  the  computation  so  that  the  mesh  overlays 
that  present  at  the  start  of  the  cycle.  A  geometric  interpretation  is  em¬ 
ployed  coupled  with  conservation  rules. 


The  first  step  is  to  decide  whether  to  add  a  new  zone.  Zones  are  added 
in  order  to  keep  the  zones  from  becoming  too  large  as  the  system  expands  due 
to  intumescence.  The  amount  of  grid  expansion  during  the  present  time  cycle 
is  simply: 


<Al)«ip»ad  ‘  "  <‘".)<dx> 


(B-37) 


So  long  as  (Ax)  .  is  less  than  Ax/2,  the  last  zone  Ax.  is  given  by: 
expand  o 


Ax  °  Ax  +  (Ax)  .  (Ax)  .  <  Ax/2 

b  expand  expand 


(B-38) 


However,  if  (Ax)  .  is  greater  than  or  equal  to  Ax/2,  then  a  new  zone  is 

expand 

added  to  the  grid: 


kmax  kmax+1  if  <Ax)expaild  1  Ax/ 2  (B-39) 

This  procedure  effectively  requires  that  the  last  zone  be  between: 


Once  the  number  of  zones  of  the  new  grid  have  been  determined,  the 
rezoning  process  begins.  Referring  to  Figure  B-7,  the  rezoning  procedure 
defines  the  variable  F^  at  x^  (e.g.,  density,  or  temperature,  etc.)  in  terms 
of  the  G^'a  at  y^  associated  with  the  expanded  grid.  From  Figure  B-7,  it  can 
be  seen  that: 


f  < 


3  Ax 
2 


(B-41) 


Expanded 


♦I«l 


L*>: 


m\ 


Coating  Did  Mot 
Expand  In  This 
Region 


Figure  B— 7 .  Schematic  of  Rezoning  Geometry 


Fk  Ax  «  Gk  Ayk  +  (cross-hatch  region)  -  (dotted  region)  (B-42) 

Fk  '  [°k  (**«  -  *t)  *  °*l  (7k  -  *k)  -  0k  (Titl  -  Vl)]/4*  <B-43> 

Expression  (B-43)  is  representative  of  the  general  procedure.  Other  cases, 
for  example  where  xfc  -  yfc  but  xfc+1  <  j  ,  are  solved  similarly. 

Care  must  be  exercised  in  rezoning  some  of  the  variables.  Where  a  grid 
point  variable  represents  the  "condition"  of  the  entire  zone,  such  as  density 
or  velocity,  the  foregoing  geometric  interpretation  is  adequate.  However,  a 
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■mss  variable  [such  as  (■  ).  ,  (m  ) .  ,  etc.]  is  computed  by  integrating  the  den- 

O  E  C  K 

sity  over  the  width  of  the  zone.  To  rezone  such  a  variable  first  requires 

transformation  to  a  "differential"  variable,  e.g.,  (■  )./(y.  -  y.).  After 

o  k  1l+i  k 

rezoning,  the  mass  associated  with  the  zone  is  found  again  by  integrating  the 

"differential"  variable  over  the  width  of  the  zone:  (m  ).  Az.  Likewise,  in- 

o  t 

stead  of  rezoning  the  temperature,  the  specific  energy,  C  T,  is  rezoned,  and 

P 

then  the  temperature  recomputed. 

B.3.9  Time  Step  (Stability  Criteria) 

The  explicit  formulation  that  has  been  developed  is  numerically  stable 
so  long  as  the  time  step.  At,  meets  certain  criteria.  A  time  step  is  computed 
to  insure  stability  of  the  energy  equation.  At  interior  grid  points,  the  con¬ 
dition  for  stability  is: 


At  <  (AX)2 
Atk  -*■  2 


Ml 


(B-44) 


where  depends  npon  the  temperature.  Thns,  the  zone  with  the  largest 
specifies  the  time  step.  The  boundary  zone,  because  of  its  different  zone 
size,  has  a  time  step  At.  associated  with  it: 


The  time  step  is  the  minimum  of  these 


At.  ,  At 
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To  catch  the  sudden  jump  at  t*0  (i.e.,  the  heat  flux  is  suddenly  "turned  on"), 
the  first  cycle  time  step  is  reduced  by  a  factor  t^  which  is  generally  of  the 
order  of  0.01: 

At0  -  tr-At  (B-48) 

The  time  step  is  then  allowed  to  get  larger  in  succeeding  cycles  but  only  at  a 

prescribed  growth  rate,  t  .  Thus,  the  time  step  for  the  computational  cycle 

S 

is  taken  as: 

Ati+1  -  tfMIN  ^t  -At1, At 1+1 1  (B-49) 


where  0  <  t^  <  1.0. 

The  continuity  equation  also  has  a  time  step  stability  criteria: 


MIN  [Ax,Ax  ] 

At  1 - “ - -  (B-50) 

C  T 

max 


where  is  the  maximum  velocity  within  the  grid  (i.e.,  the  velocity  of  the 

boundary,  tmt  **  vVmjt)  .  After  determining  the  velocities  and  densities,  the 
program  checks  for  consistency  of  the  current  time  step  [Equation  (B-46)]  to 
insure : 


At  1  At  (B-51) 

6 

If  this  relationship  is  not  satisfied,  the  program  computes  a  new  At  ■  t^-Ato 
and  then  recycles  through  the  energy  equation,  etc. 
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B.3.10  Tabular  Entries  and  Search  Subroutines 


The  data  as  a  function  of  temperature  for  enthalpy,  specific  heat,  and 
chemical  heat  production  have  all  been  entered  into  the  computer  program  as  a 
subroutine  called  BLOCK  DATA.  The  values  for  each  of  the  individual  constitu¬ 
ents,  Part  A,  Part  B,  and  Part  A+B,  have  been  entered.  Then  each  of  the  three 
thermodynamic  functions  are  found  for  the  appropriate  temperature  by  a  search 
routine  which  linearly  interpolates  between  values  (where  applicable),  i.e.,: 


(B-52) 


where  6^  and  are  thermodynamic  variables  at  the  respective  tempera¬ 

tures  Tj  and  Ti+^,  and  F  is.  the  interpolated  value  (0,^  <.  F  <  9i+1)  for  the 
temperature  T. 

If  the  coating  is  assumed  to  be  a  composite  of  the  sum  of  its  constitu¬ 
ents,  then  the  thermodynamic  quantity  is  F,  computed  as  a  sum  of  the  individ¬ 
ual  constituent's  thermodynamic  quantity,  Fj,  times  its  respective  concentra¬ 
tion,  Cj,  in  the  mixture: 


X>i'j 

j-1 


(B-53) 


where  there  are  J  total  constituents.  The  Cj's  *  function  of  temperature 
and  are  determined  by  using  Equation  (5-22). 


B.3.11  Thermal  Conductivity 


There  is  a  thermal  conductivity  associated  with  the  virgin  material, 

K  ,  and  another  values  associated  with  the  char,  K  .  The  program  provides  a 
o  c 

smooth  transition  between  these  two  values  according  to  the  density  change: 


i 


i 


B-25 


DISTRIBUTION  LIST  (Continued) 
REPORT  NO.  NADC-822 11-60 


U.  S.  Polymeric,  Inc.  (D.  Beckley) ,  700  East  Dyer  Road,  Santa 

Ana,  CA  92707  . 

University  of  Dayton  Research  Institute  (I.  Salyer),  300  College 
Park  Avenue,  Dayton,  OH  45469  . 

Wright  Aeronautical  Laboratories  (MLBE,  MLSA) ,  Wright-Patterson 
AFB,  OH  45433  . 


